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^ ■ We present a general all-order formulation of Sudakov resummation in QCD in terms 
of dispersion integrals. We show that the Sudakov exponent can be written as a dis- 
persion integral over spectral density functions, weighted by characteristic functions that 
encode information on power corrections. The characteristic functions are defined and 
computed analytically in the large-/?o limit. The spectral density functions encapsulate 
the non-Abelian nature of the interaction. They are defined by the time-like discontinuity 
of specific effective charges (couplings) that are directly related to the familiar Sudakov 
anomalous dimensions and can be computed order-by-order in perturbation theory. The 
dispersive approach provides a realization of Dressed Gluon Exponentiation, where Su- 
dakov resummation is enhanced by an internal resummation of running-coupling correc- 
tions. We establish all-order relations between the scheme-invariant Borel formulation 
and the dispersive one, and address the difference in the treatment of power corrections. 
We find that in the context of Sudakov resummation the infrared-finite-coupling hypoth- 
esis is of special interest because the relevant coupling can be uniquely identified to any 
order, and may have an infrared fixed point already at the perturbative level. We prove 
that this infrared limit is universal: it is determined by the cusp anomalous dimension. 
To illustrate the formalism we discuss a few examples including B-meson decay spectra, 
deep inelastic structure functions and Drell-Yan or Higgs production. 
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1 Introduction 



An accurate theoretical description of inclusive differential cross sections and decay spec- 
tra is essential for many aspects of collider and flavor physics; it is therefore one of the 
primary goals of QCD perturbation theory. The main challenge in describing differential 
distributions is associated with kinematic regions where there is a large hierarchy of scales, 
as occurs for example in the production of a jet with a large energy and a small mass. Near 
the exclusive limit of phase space, the so called Sudakov region, only soft and collinear 
radiation is kinematically allowed, and because of the dynamical enhancement of such 
radiation, multiple gluon emission becomes important. In this way the hierarchy of scales 
translates into large logarithms that appear with increasingly high powers in perturbation 
theory and spoil the convergence of the expansion. Even a qualitative description of the 
distribution in this region require all-order resummation. Classical examples where the 
Sudakov region has been thoroughly studied include deep inelastic structure functions for 
x — > 1 [1-5], Drell-Yan and Higgs production near the partonic threshold [1,2,6-13] or at 
small transverse momentum (see e.g. [14-18]), event-shape distributions [19-29], heavy- 
quark fragmentation [30-34], and inclusive B decays into a light quarks [32,35-44], e.g. 
B -> X s 7 or B -> XJv. 

From a theoretical perspective the Sudakov limit is very interesting. Despite the fact 
that one is dealing with a complex multi-scale problem in a non-Abelian gauge theory, 
one can in fact control the dominant perturbative corrections to all orders and resum the 
series. The simplification of the expansion in the Sudakov limit is a direct reflection of 
the factorization property of infrared (soft and collinear) singularities. In infrared safe 
quantities the singularity associated with soft and collinear radiation (integrated over 
phase space) cancel exactly, order by order in a s , with infrared singularities in virtual 
corrections. The structure of these singularities and their cancellation is encoded in the 
resummation formalism: the large logarithms are nothing but the finite remainder in the 
sum of the real and virtual contributions, which are separately divergent. 

Sudakov resummation is best formulated in moment (Mellin) space where the multi- 
gluon phase space factorizes. In moment space the combined effect of soft and collinear 
radiation and the corresponding virtual corrections appears an exponential Sudakov factor. 
The Sudakov factor sums up to all orders the dominant radiative corrections that are 
enhanced by powers of In N (where N is the moment index), while neglecting corrections 
that are suppressed by powers of N for iV — > oo. 

While Sudakov resummation is designed to deal with logarithmic singularities alone, 
infrared sensitivity appears also through power-suppressed effects, O ((A/m) p ) [22,45-53], 
where m is the hard scale and where p is an integer. In the Sudakov region, power 
corrections are particularly important: whereas perturbative corrections are enhanced at 
iV — > oo by powers of the logarithm In TV, power corrections are enhanced by powers 
of N, taking the form O ((NK/m) p ). This parametric enhancement of non-perturbative 
corrections implies that for any fixed coupling (or a fixed hard scale) the N — > oo limit 
itself is strictly beyond the reach of resummed perturbation theory. Put differently, upon 
ignoring these power corrections one can only expect the large-iV resummation to be a 
valid approximation in a restricted range: N <C m/A. 
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Quite remarkably, perturbation theory itself is sensitive [47] to the presence of these 
power-suppressed corrections [54] . This sensitivity appears through infrared renormalons 
[55,56], factorially increasing coefficients that emerge out of the integral over the running 
coupling. Despite the fact that the Sudakov exponent can be uniquely computed to any 
logarithmic accuracy, the series as a whole diverges. The understanding that renormalons 
are an inherent part of the Sudakov exponent [47,52,53] has prompted the development of 
a new resummation formalism, Dressed Gluon Exponentiation [22,23,57]; for a recent 
review, see [54]. DGE combines Sudakov resummation with an additional "internal" 
resummation of running-coupling corrections^. The internal resummation is based on 
an all-order calculation of the exponentiation kernel in the large-/3o limit, which exposes 
the renormalons and thus opens the way for including non-perturbative power corrections. 

DGE has been applied and successfully compared with data in a variety of inclusive dis- 
tributions: event-shape distributions [22,23], deep inelastic structure functions [5], heavy 
quark fragmentation [30] and inclusive decay spectra [40-44]. It has proven to provide a 
good description of these distributions over the entire Sudakov region, significantly extend- 
ing the range of applicability of resummed perturbation theory. In its first application, to 
event-shape distribution, DGE has opened the way for a quantitative description of the 
two-jet region, which had been inaccessible to analytic methods before. In heavy-quark 
fragmentation and inclusive decay spectra DGE provides a viable alternative to the con- 
ventional approach where the corresponding distributions ("heavy-quark fragmentation 
function" and "shape function", respectively) have been parametrized by some ad hoc 
functional forms. Further details can be found in recent review talks [44, 54] and in the 
original publications. 

Technically these results have been obtained by trading [22] the integration over the 
running coupling in the Sudakov exponent for a scheme-invariant [64] Borel integration. 
The Borel transform has been computed analytically in the large-/3o limit, allowing to 
identify the renormalon singularities. The perturbative sum is then defined by the Prin- 
cipal Value prescription and power corrections are parametrized based on the renormalon 
ambiguities. While the original calculations of event-shape distributions in Refs. [22, 23] 
have been performed at next-to-leading logarithmic (NLL) accurac}{§ it was observed that 
the Borel formulation is in fact completely general, and can accommodate subleading cor- 
rections with any logarithmic accuracy, provided that the relevant Sudakov anomalous 
dimensions are known. This has been consequently used in other applications [5,32,43,54] 
to perform calculations with formal NNLL accuracy while preserving the pattern of renor- 
malon singularities found in the large-/?o limit. 

In this paper we establish an alternative formulation of Sudakov resummation that 

* "Internal" refers here to the fact that the argument of the relevant logarithms involves loop-momenta 
that are integrated over, rather than external scales. Resummation of running-coupling corrections has 
been extensively studied in the past, primarily in the context of single-scale quantities [45,49,55,56,58-62]. 
See in particular the discussion in section 7 of [58] and in [45] . For a related early discussion of Sudakov 
resummation for the electron form factor in QED see also [63]. 

§A proper definition of the QCD scale A corresponding to the "gluon bremsstrahlung" coupling [12] 
(see Eq. (104) in [22]) is sufficient to account for the entire non-Abelian correction, 0(Ca/ fto), at this 
order; A is fixed by the NLO correction to the cusp anomalous dimension [65,66]. 



3 



captures the same class of radiative corrections as well as the renormalon structure, using 
the dispersive approach^. Following [68, 69] we show that the Sudakov exponent can be 
written, to any logarithmic accuracy, as a dispersive integral of the time-like discontinuity 
of specific non-Abelian effective charges, integrated with respect to characteristic functions 
that are defined (and computed analytically) in the large-/?o limit. This way the exponent 
in the non-Abelian theory is written in full analogy with the way renormalon resummation 
has been formulated in Refs. [49,61,62] at the level of a single dressed gluon (the large-/? 
limit) with one important difference: the "effective" coupling is computed order-by-order 
in the non-Abelian theory. 

Refs. [68,69] emphasized the fact that Sudakov resummation leaves some freedom (see 
also [70]) that does not allow to make conclusive statements about power corrections. 
Indeed, in general there is a conflict between the heuristic argument for power corrections 
in expressions involving integration over the running coupling (e.g. [47]) and renormalon- 
based arguments [52]. It was observed that, at the perturbative level, the resummation of 
Sudakov logarithms can be implemented in infinitely many ways, each of which consisting 
of an integral mapping of a perturbative object (the 'running Sudakov coupling') onto a 
fixed target: the 'Sudakov exponent'. To each mapping corresponds a different kernel. In 
Ref. [68] it was proposed to fix this ambiguity and select the unique resummation scheme 
by requiring that in the large-/9o limit, the 'running Sudakov coupling' coincides with a 
dressed gluon propagator. The specified scheme so determined is the dispersive approach, 
where the above conflict does not arise. In the present paper we further develop this 
approach, establishing the connection with the physical Sudakov anomalous dimensions. 

We show that this dispersive formulation is in close correspondence with the scheme- 
invariant Borel formulation and establish the relations between the two. We further show 
that the effective charges appearing in the dispersive formulation are uniquely defined to 
any order in perturbation theory and stand in one-to-one correspondence with the relevant 
Sudakov anomalous dimensions. These effective charges provide a (process-dependent) 
generalization of the concept of the "gluon bremsstrahlung" coupling, which was defined 
in [12] at the next-to-leading order. Moreover, we find that the dispersive formulation is 
closely related to the joint resummation formalism of Ref. [17], where the effective charges 
are given a direct diagrammatic interpretation in terms of "webs" [71]. 

The dispersive approach proposed here is particularly suited to the implementation of 
the infrared finite coupling hypothesis of Ref. [49]. This possibility is of special interest 
in the context of Sudakov resummation, since the relevant non-Abelian coupling can be 
identified to any order, and furthermore, may have an infrared fixed-point already at the 
perturbative level. In this case the coupling may have a causal analyticity structure [72-74] , 
free of any Landau singularities. 

We proceed as follows: in Sec. [2] we recall some general features of the Sudakov limit 
and define the basic elements used in Sudakov resummation. Then, in Sec. [3] we analyze 
the structure of the exponentiation kernel in momentum space. We first use general 
considerations of infrared factorization and renormalizat ion-group invariance to identify 

dispersive representation of the Sudakov exponent was first considered in the original DGE pa- 
per [22], see Eq. (43) there. In that paper, however, the discussion was limited to NLL accuracy. 
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the physical Sudakov anomalous dimensions. We then employ all-order resummation in 
the large-/?o limit arriving at a dispersive representation of the kernel. Next, in Sec. H] 
we generalize the dispersive representation to the full non-Abelian theory. We show that 
this leads to a unique (yet process-dependent) generalization of the concept of a "gluon 
bremsstrahlung" coupling, the "Sudakov effective charges", which are computed order- 
by-order in perturbation theory. We then analyze the evolution of these effective charges, 
finding that they may have a perturbative infrared fixed point. We further prove that 
the corresponding fixed-point value is process-independent, and compute its Banks-Zaks 
expansion. Then in Sec. [5] we use the tools of the previous sections to derive a general 
dispersive representation of the Sudakov exponent. We also present explicit results for the 
characteristic functions in several different processes and analyze their properties. In Sec. [6] 
we explain the relation between the dispersive formulation and the scheme invariant Borel 
formulation. Next, in Sec.[7]we address power corrections. Finally, in Sec.[8]we summarize 
our conclusions. 



General set-up: Sudakov resummation for inclusive 
distributions 



Consider an inclusive infrared-safe distribution, dT(m 2 ,r)/dr where the r 
characterized by Sudakov logarithms, 

dr(m 2 , r) 



limit is 



dr 



6(r)[l + 0(a 3 )) + C 



7T 



ln(r) b\ — di 



+ regular terms > + 



(2-1) 

whose maximal power grows with order as a s n In 2 " 1 (r)/r due to multiple soft and collinear 
radiation. The notation hi and di used for the coefficients is associated with their phase- 
space origin; this will be elucidated in the next section. The integration prescription [] + 
is defined bjU 

1 1 £ ^ (F(r) - F(0)) , (2.2) 



drF(r) 



where F(r) is a smooth test function. This prescription accounts for the divergent virtual 
corrections, which cancel against the singularity generated when integrating the real- 
emission contributions near r = 0. 

In order to describe the distribution at small r we need to resum all the singular terms 
in (12.11) to any order. To this end it is convenient to work in moment space, where the 
multi-gluon phase space factorizes. The moments are given by: 



T N (m 



dr(l 



_ : dT(m 2 ,r) 
dr 



(2.3) 



where the plus prescription in (12.11) guarantees that the moments are finite, cf. (\2.7\i 
below. Once resummation has been performed analytically in moment space one recovers 



II Differentiation with respect to u can be used to generate powers of ln(r). See Eq. (6.4) in [75] for a 
more general definition. 
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the distribution in momentum space by an inverse Mellin transformation: 



dY{m 2 



dr 



rcsummed 



loc dN 



where the integration contour runs parallel to the imaginary axis in the complex N plane, 
to the right of the singularities of the integrand. 

We assume that the large- N limit singles out infrared singularities associated with two 
physical scales, which we call the jet scale — momenta of order 0(m 2 r) — and the soft 
scale — momenta of order 0(m 2 r 2 ). We therefore have a double hierarchy of scales: 

m 2 > m 2 r > m 2 r 2 < > m 2 > m 2 /N > m 2 /N 2 . (2.5) 

While we are interested in large N, in order to apply perturbation theory we assume that 
all three scales are within the perturbative regime, namely m 2 /N 2 ^> A 2 . 

This general scenario, with minor variations, is encountered in many different applica- 
tions, for example: 



Inclusive B decays into light quarks [32,35-44], e.g. B — > X s ^ oqI? — ► X U W. These 
decays are characterized by jet-like kinematics having a large hierarchy between the 
energy of the hadronic system X and its mass, or equivalently, between the two 
lightcone components of the jet, p ± = E =p namely p~ ^> p + . Thus, in (12. 5p 
the hard scale is the b-quark mass or the large lightcone component of the jet 
(p~), while r is the ratio between the small and large components r = p + /p~ . The 
intermediate scale is the jet mass p + p~ ~ m 2 r. The 'soft' scale m 2 r 2 is associated 
with soft gluons that couple of the b quark, putting it slightly off its mass shell 
prior to the decay. In perturbation theory the soft subprocess is therefore associated 
with the momentum distribution of the b quark inside an initial on-shell b quark, 
see [32] and Refs. therein. Beyond perturbation theory this function is replaced 
by the momentum distribution of the b quark in a meson, which differs from its 
perturbative counterpart by power corrections. 

Inclusive B production in e + e~ annihilation, see e.g. [30-34]. Here r = 1 — x where 
x = 2Eb/Q, the energy fraction of the detected B meson and Q is the center-of-mass 
energy. At a difference with (I2.5P the three scales are Q 2 3> Q 2 {1 — x) ^> m 2 (l — x) 2 , 
where m& is the b-quark mass. The soft scale m 2 r 2 is associated with the heavy quark 
fragmentation process (which, at the perturbative level, is similar to b decay [32]). 

Drell-Yan or Higgs production in hadronic collisions in the DIS scheme [1,2,6-12], 
where m = Q is the mass of the produced pair (or ran, in case of Higgs production), 
r = 1 — Q 2 /s is the fraction of energy carried by soft gluon radiation into the final 
state, and s is the partonic center-of-mass energy. Here the soft scale is genuinely 
associated with the Drell-Yan process while the jet mass scale enters via the quark 



**Note that in case of the semileptonic decay b — > uW, resummation is applied to the triple differential 
width [42] . For simplicity of the notation we do not write here explicitly derivatives with respect to other 
kinematic variables. 
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distribution function in the DIS factorization scheme (the x — > 1 limit of deep inelas- 
tic structure functions) and it can be traded for \j? F upon using the MS factorization 
scheme. 

• Event -shape distributions in e + e~ annihilation, near the two-jet limit [19-25,29]. 
Here r is identified with a shape variable that vanishes in the two-jet limit, for 
example 1 — T where T is the thrust variable or C/6 where C is the "C parameter" ; 
the hard scale m is the center-of-mass energy in the collision. Because jet observables 
are less inclusive [22,23,29,76], the formalism we develop applies only to a limited 
logarithmic accuracy. Non-inclusive effects can be accommodated in a class of event- 
shape variables (where Eqs. (12. 7p and (12.81) below are valid [19]) but this goes beyond 
the scope of the present work. 

As we shall see the same resummation formalism applies in all these examples. Moreover, 
the coefficients of the logarithms have a certain degree of universality: first, the leading 
logarithms are always related to the cusp anomalous dimension [65,66], A(a s ) in (12. 8p 
below; second, the physics on the jet-mass scale m 2 r is the same in all these examples: 
it is associated with the production of a quark jet with a constrained invariant mass of 
0(m 2 r). The corresponding Sudakov logarithms are the same as in deep inelastic structure 
functions in the x — > 1 limit; they are generated, in all these examples, by the anomalous 
dimension function B(a s ) in (12. 8p . In contrast, the physics on the soft scale m 2 r 2 is rather 
different; this will be reflected in the resummation formalism through a process-dependent 
Sudakov anomalous dimension V(a s ). 

Owing to the factorization property of soft and collinear radiation the moments (12. 3p 
can be written as follows: 



where the Sudakov factor Sud(m 2 , N) sums up to all orders in perturbation theory all the 
terms that diverge for N — > oo, H(m 2 ) is a hard coefficient function that includes constant 
(TV-independent) terms at each order in a s ; such corrections arise from virtual diagrams, 
proportional to S(r). Finally, AT^lm 2 ) = 0(1/N) includes residual real-emission correc- 
tions that fall at large N as 1/N (up to logarithms). In this paper we are interested in the 
Sudakov factor Sud(m 2 , N). Both H(m 2 ) and AT^(m 2 ) can be computed order by order 
in a s by "matching" the resummation formula with the fixed-order expansion. We shall 
not discuss them further here. 

The general formula [1,2] (see also [3,19,31,32,38,41-43]) for the Sudakov factor is: 



where the subtraction of 1 accounts for the virtual corrections encapsulated in the plus 
prescription in (I2.ip . and the momentum-space kernel, which is fully defined by the real- 
emission contributions that are singular for r — > 0, takes the general form: 



T N (m 2 ) = H(m 2 ) x Sud(m 2 , N) + AT N (m 2 ) , 



(2.6) 




(2.7) 



C R 



R(m 2 ,r) 1 



r 



r 




2 



2 dtf 

2 k 2 



A{a s {k 2 )) + B (a s (rm 2 



)) - V (a s (r 2 m 2 



)) 



(2.8) 
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where A, B and T> are Sudakov anomalous dimensions that can be computed order-by- 
order in a s . A(a s ) is the universal cusp anomalous dimension [65,66], whose expansion, 
in the MS scheme is known to three-loop order [77]: 



■*(<».) = c f 

Po 



a + a 2 a 2 + a 3 a 3 + 



(2.9) 



where a = [3o(xf s /n. The other Sudakov anomalous dimensions appearing in Eq. (I2.8P are 

Po 

which is associated the the jet mass scale m 2 r, and 



bid + b 2 a 2 + 6 3 a 3 + ■ ■ 



(2.10) 



Po 



^3 



(2.11) 



which is associated with the soft scale m 2 r 2 . The known expansion coefficients of these 
anomalous dimensions are summarized in Appendix |A] As already mentioned, T> is process 
dependent. In Appendix |A] we quote the coefficients in two examples, one of Drell-Yan (or 
Higgs) production and the second of inclusive B decay (or heavy-quark fragmentation - 
the anomalous dimension is the same [32]). These two examples will be used throughout 
this paper. 

Next, note that Cr, the overall coefficients of all the anomalous dimensions entering 
(12. 8p in a given process, depends on the color charge of the hard parton(s) that radiates. In 
inclusive B decay or heavy-quark fragmentation Cr = Cf, while for Drell-Yan production 
Cr = —2Cf and for Higgs production by gluon-gluon fusion Cr = —2Ca- 

Finally, it is useful to introduce some alternative definitions for the Sudakov factor 
using the Laplace weight instead of the Mellin one. The first observation is that at large 
iV one can make the following approximation of the real-emission weight factor in (12.71) 

(1 _ r f-l „ e -Nr ? 

which only modifies the exponent by 0(1/ N) terms. Next one notes that as far as the 
real-emission part is concerned the integral over r can be extended to r = oo, well beyond 
the physical phase space r = 1: the resulting modification of the Sudakov exponent is 
exponentially small, 0(e~ N ). This, however, does not apply to the virtual terms, where 
any change of the upper limit of integration is reflected in a modification of the constant 
O(N ) term for iV — ► oo. Let us therefore define: 

dr 



Sud(m 2 , N) = exp \C R — 



-Nr 



9(r < 1) R(m 2 



which satisfies 



In ( Sud(m 2 , N)j = In (Sud(m 2 , N)) + 0(l/N). 



(2.12) 



(2.13) 



The 0(1/ N) contributions by which Eq. (12.121) differs from Eq. (12.71) . can be compensated 
in the expression for the moments by a different remainder function. The physical moments 
T^(m 2 ) can be obtained by "matching" with the fixed-order result: 



T N (m 2 ) = H(m 2 ) x Sud(m 2 , N) + AT N (m 2 



(2.14) 
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where H(m 2 ) is the same hard function as in (12.61) . while Ar^^ 2 ) can be determined to 
fixed order from the perturbative expansion of T^m 2 ). 

It is useful to further define a Sudakov factor where the subtraction term too is inte- 
grated to r = oo: 



The suppression of the weight factor e~ Nr — e~ r guarantees convergence for r — > oo. Again, 
the manipulations used in obtaining (12.151) from (12. 7p do not change the N —>■ oo divergent 
terms, which all emerge from the r — > limit. Here however, in contrast with (I2.12p . the 
constant terms O(N ) for N — > oo are modified. 

Note that the moment-space exponent of Eq. (12. 7p is composed of harmonic sums 
(powers of the \&(iV) function, which differ from ln(A^) by constants and inverse powers 
of N) whereas the Laplace version, Eq. (I2.15p . involves strictly powers of ln(iV), with no 
additional constants nor 0(1/N) terms. Eq. (I2.12p contains powers of ln(iV) and constant 
terms, with no additional 0{l/N) terms [68,78]. Both fl2~T51 and (1277]) . but not fl2~T2l . 
are defined such that for N — 1 the exponent exactly vanishes so the Sudakov factor 
becomes 1. 

Both the constant term, O(N ), and the 0(1/N) contributions by which Eq. (I2.15P 
differs from Eq. (12.71) . can be compensated in the expression for the moments by a different 
multiplicative hard function and a different remainder function: 



similarly to H(m 2 ) and ATn{itl 2 ) in (12. 6ft . H(m 2 ) and AT^im 2 ) in (12.161) can be de- 
termined to fixed order from the perturbative expansion of rjv(m 2 ). Finally, Eq. (12.161) . 
similarly to (12.141) and (12.61) . can be inserted into the inverse Mellin transformation (12 .4p 
to generate the resummed distribution in momentum space. 

3 The exponentiation kernel 

The key to constructing an effective resummation technique is the understanding of the 
all-order structure of the kernel. The purpose of this section is to analyze the momentum- 
space kernel (12.81) . first using the general considerations of infrared factorization and 
renormalizat ion-group invariance, and then by considering the all-order calculation of 
the kernel in the large-/5o limit. 

3.1 Infrared factorization and renormalization— group invariance 

Let us begin by noting that Sudakov logarithms in infrared and collinear safe distributions 
arise from different regions of phase space, involving different physical scales. As follows 
for (12.11) . in the perturbative expansion of the kernel, 




(2-15) 



T N {m 2 ) = H{m 2 ) x Sud(m 2 , N) + AT N (m 2 ); 



(2.16) 



C R 



R(m 2 , r) 



C R 



ln(r) bi — di a s (m) 



+ ... 



(3.1) 



r 
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the two are mixed together. In contrast, upon taking one derivative (cf. Eq. (I2.8p ) 



dR [ m2 ^ = J(rm 2 )-S(r 2 m 2 ), (3.2) 
dm m z 

the dependence upon the 'jet' (rm 2 ) and 'soft' (r 2 m 2 ) scales is nicely disentangled into 
two different functions corresponding to the two subprocess, J and S respectively. These 
are the physical 'jet' and 'soft' Sudakov anomalous dimensions. They are related to the 
ones appearing in (12.81) by the following relations: 

C R J{^)=A{a s {^)) + dB ^ 2 2 )) 

a in /i z 

C R S(S)=A M + d ^0». (3 ' 3> 

We emphasize that j7(/i 2 ) and S(fi 2 ) are renormalizat ion-group invariant functions, while 
their conventional decomposition into the cusp A, B and T> (which are renormalizat ion- 
scheme-dependent quantities) is specific to the MS scheme. 



3.2 The physical interpretation of the Sudakov anomalous di- 
mensions 

The functions J{[i 2 ) and S(fi 2 ) have a clear physical meaning: they govern the scale- 
dependence of certain physical quantities (the 'physical evolution kernels' or 'physical 
anomalous dimensions') in the Sudakov limit. Let us illustrate this statement using a 
couple of examples. 

Deep inelastic structure functions 

The scale-dependence of the (flavour non-singlet) deep inelastic structure function F 2 
can be expressed in terms of F 2 itself, yielding the following evolution equation (see 
e.g. Refs. [63,79,80]): 



dF 2 (x,Q 2 ) f 1 dz 
dlnQ 2 



f dz 

/ -K ms(F2) (x/z,Q 2 )F 2 (z,Q 2 ). 

J X z 



(3.4) 



Kms(F 2 ){ x i Q 2 ) is the physical evolution kernel; it is renormalization-group invariant. 
Defining moments by 

F 2 (N,Q 2 )= [ dxx N ~ l F 2 {x ) Q 2 )/x ) (3.5) 
Jo 

Eq. (13. 4p implies that the moment-space physical evolution kernel, the 'physical anomalous 
dimension', is: 

K ms{F2) (N,Q 2 ) = J" dxx N - l K mB{F2) {x ) Q 2 )/x = ^J^g/^ - (3-6) 
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Let us now consider the x — > 1 limit. In this limit the evolution of the structure 
function takes a simple from [4,5] (see Eq. (13) in [54]): 



d\nF 2 (N, Q 2 
d\nQ 2 



„N-1 



dx- 



1 



x 



j((l-x)Q 2 ) + 



dlnH(Q; fMp) 



+ 0(1/N) (3.7) 



d\nQ 2 

where the first term, which includes the N —>■ oo divergent corrections to all orders, is 
controlled by the Sudakov anomalous dimension J7"(/i 2 ) of Eq. (I3.3p . The constant term 
can be written [70] in terms of the quark electromagnetic form factor ^(Q 2 ) [81-84]: 



d\nH(Q; fi F ) d\n(T(Q 2 )Y 



dlnQ 2 



rflnQ 2 



+ C f 



dp 
~7 



J(S 



(3.f 



G(l,a s (Q 2 ),e = 0) + B(a s (Q 2 )) , 



where each of the two terms in the first line is separately infrared divergent, but the 
divergence cancels in the sum [70] ; in the second line the result is expressed in terms of B of 
Eq. (12.1 Op above and G (Q 2 /fi 2 , a s (// 2 ), e), which is the finite part of <iln (F(Q 2 )) 2 / dhiQ 2 
as defined in Ref. [84] using dimensional regularization. 

Comparing (13. 6p and (13.71) we therefore find the following relation in momentum space: 



DIS(F 2 ) 



>,Q 2 



a 



J((l-x)Q 2 
1 — x 



J{{l-x)Q 2 
1 — x 



d\n(F(Q 2 )) 2 
dlnQ 2 

( 



5(1 - x) + O ((1 - x)°) 



dln(F(Q 2 )) 2 
dlnQ 2 



+ Ci 



Q 



dp 



infrared finite 



6(1 



x 



+ 0((l-x) ) , 

(3.9) 

where the plus prescription is defined in ( 12.21) and the O ((1 — x)°) term neglected here 
is integrable for x — > 1 (although divergent: it contains powers of ln(l — x)). We thus 
find that Cp J ({1 — x)Q 2 )/(l — x) is the leading term in the expansion of the physical 
evolution kernel K D1B tp 2 -j(x, Q 2 ) in the x — > 1 limit with fixed jet mass (1 — x)Q 2 . The next 
term in this expansion, proportional to 5(1 — x), is comprised of purely virtual corrections 
associated with the quark form factor. This term is infrared divergent, but as indicated in 
the second line in ( I3.9p . the singularity cancels exactly upon integrating over x with the 
divergence of the integral of CpJ ((1 — x)Q 2 )/(l — 

Drell-Yan cross section 

The cross section of the Drell-Yan process, h a + 



x) near x 



1. 



e 1 e 



X, is: 



da 



Ana' 



dQ 2 9Q 2 s 



E 



1,3 



dx i dx •) 



X 3 



fi/h a (Xi,^ F ) f. 



j/h h { x ji 



Hf) gij(r,Q 2 ,fJ.p) 



(3.10) 



where s is the hadronic center-of-mass energy, s = XiXjS is the partonic one, Q 2 is 
the squared mass of the lepton pair and r = Q 2 /s. The partonic threshold r — ► 1 is 
characterized by Sudakov logarithms. 
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Let us define the Mellin transform of the quark-ant iquark partonic cross section in 
(I3.10p . e.g. in electromagnetic annihilation g qq (r, Q 2 , fi 2 F ) = e 2 [5(1 — r) + 0(a s )], by 



G qq (N, Q 2 , n*) ee f dr r N ~ l g qq (r, Q 2 , /4) . 
Jo 

Considering the Sudakov limit one has 

G qq (N, Q 2 , nl) ~ H DY (Q 2 , n%) Sud DY (iV, Q 2 , + 0(1/N) 

with 



(3.11) 



(3.12) 



Sud DY = exp < 2 / dr 



1 ^N-l 



1 - T 



dk 



I-'f 



A (a s (k 2 )) ^— + X? DY («.((! - rfQ 2 )) 



(3.13) 



Upon taking the logarithmic derivative we identify the physical anomalous dimension iS DY 
of Eq. (133|) : 



din G qq (N,Q 2 ,/j 2 F ) 
dlnQ 2 



K DY (N,Q 2 ) = I drr N - l K m {T,Q 



2C F / dr 



i t n-i _ ! 



1 -r 



5 DY (Q 2 (l - r) 2 ) + + {l/N) , 



dlnQ 2 



(3.14) 



where K uy (t, Q 2 ) is the physical Drell-Yan evolution kernel defined for arbitrary r; in the 
second line we considered the large- N limit where Eqs. (I3.12p and (I3.13P apply and then 
used Eq. ( jSP - 

The constant term can be expressed in terms of the analytically-continued electromag- 
netic quark form factor [70] (see also [6,8,9]): 



dln# DY (Q 2 ,/4) dln\(F(-Q 2 ))\ 2 



dlnQ 2 



dlnQ 2 



+ C t 



Q 



o F 



S m {n 2 ) 



(3.15) 



where the infrared singularities cancel in the sum, as in (13.81) . Thus, in momentum space 
we have: 



2\\|2 



K DY (r,Q 2 )=2C F 



+ 



2C F 
( 



S DY (Q 2 (l-r) 2 ) , dln\(F(-Q 2 ))\ 

1 - T 

S» y (Q 2 (1-t) 2 



dlnQ 2 



5(1 - t) + O ((1 - r)°) 



+ 



dln|(^(-Q 2 ))| 



1 - r 



2\\\ 2 



(3.16) 



\ 



dlnQ 2 



d[if 
o r 



5 DY (yU 2 



v 

infrared finite 



*(1 - T ) + ((1 - r)°) 



12 



We see that the Sudakov anomalous dimension S DY controls the leading term in the expan- 
sion of the physical Drell-Yan kernel (13.141) near threshold. The r — > 1 limit is taken such 
that Q(l — t), corresponding to the total energy carried by soft gluons to the final state, is 
kept fixed. The next term in this expansion, proportional to 8(1 — t), is determined by the 
quark form factor, analytically-continued to the time-like axis. This purely virtual term 
is infrared singular, but upon performing an integral over r this singularity cancels with 
the one generated by integrating the real-emission term 2CpS DY (Q 2 (l — t) 2 ) / (l — t) near 
t -> 1. 



3.3 The large— (3q limit: a dispersive representation of the kernel 

Upon disentangling the 'hard', 'jet' and 'soft' scales in (13.21) we have significantly reduced 
the complexity of the problem: we are now dealing with two independent single-scale 
quantities, J(rm 2 ) and S(r 2 m 2 ). At this point we need additional tools to examine 
the all-order structure of these two anomalous dimensions. As these objects define the 
exponentiation kernel, an approximation based on a single dressed gluon appears most 
natural [22]; multiple emission is taken into account by exponentiation. In this section we 
shall therefore consider the single-dressed-gluon approximation, the large-/?o limit, post- 
poning the discussion of multiple emission to Sec. [5]where we shall perform exponentiation 
in moment space. 

A convenient way to compute the Sudakov anomalous dimensions in the single-dressed- 
gluon approximation is the dispersive techniqud^- A brief description of the method is given 
in Appendix [B] In order to compute the real-emission contribution to some physical quan- 
tity to all orders within this approximation one first evaluates the single gluon emission 
diagrams with an off-shell gluon, k 2 = /i 2 ^ 0. The result of this 0(a s ) calculation de- 
fines the "characteristic function" J-"(fi 2 /m 2 ,r). This function is then integrated with the 
discontinuity of the coupling, generating the all-order sum. Here we are interested specifi- 
cally in the singular contributions for r — > 0. In taking this limit we expect to identify the 
dependence on the two scale rm 2 and r 2 m 2 corresponding to the anomalous dimensions 
J and S, respectively. 

Dispersive representation of the kernel 

The result for the characteristic function at small r takes the form: 



C R ^f(e,r) = C R ^ - 

7T 71 r 



?J (e/r) - f s (e/r 2 ) 



regular terms , (3-17) 



V 



^The dispersive technique is a general method to perform resummation of running-coupling correc- 
tions. This technique was developed in Refs. [49,61,62], with no specific consideration of the Sudakov 
limit; it was used to study power corrections in a variety of applications. A well-known example is 
the average thrust [45,48,49]. In Ref. [22] the dispersive technique was used to compute the Sudakov 
exponentiation kernel for the thrust distribution based on the results of Ref. [46] for the characteristic 
function. Similar calculations using the dispersive technique were later done for other event-shape vari- 
ables including the heavy-jet mass [23], the C parameter [24] and angularities [25], as well as for heavy 
quark fragmentation [30]. 
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where e = fi 2 /m 2 , the ratio between the gluon virtuality fi 2 and the hard scale in the pro- 
cess, m 2 . The singular terms are distinguished by considering the limits lim r ^ ( r ^ r ( e ; r )) 
with either fixed e/r or fixed e/r 2 . The terms that are regular in these two limits are 
irrelevant for Sudakov resummation. A detailed example of how J- aia& .{e, r) is constructed 
by considering these two limits can be found in Sec. 3 of Ref. [24], where the kernel of the 
C parameter was computed (note the difference in complexity(l) between the full !F(e, r) 
and Fsing.ie, r)). Below we shall summarize the results for J r sing .(e, r) in some examples. 

The central point is that Eq. (13.171) exhibits the general property of infrared factor- 
ization summarized by Eq. (13. 2j) : jF(e, r), which is a function of two variables can be 
decomposed in this limit into a sum of two functions of a single variable: one depends on 
the jet mass scale and the other on the soft scale. Such a separation of scales cannot be 
done in an ordinary fixed-order calculation: it requires keeping an infrared regulatoi@ in 
place. 

With a single gluon emission the contribution to the physical quantity itself and to 
the exponentiation kernel are the same (see the arguments leading to Eq. (13.381) below). 
We can therefore use (I3.17P to construct the resummed kernel in the large-/3o limit as a 
dispersive integral according to Eq. (IB. 101) : 



C 



n 



R(m 2 ,r)_ = Cr r°°d^ 

largo fa 00 JO A*- 2 



■py(/i 



2, 



Fsinsilf/m 2 , r) - .F sing .(0, r) 



(3.18) 



where we neglected the power corrections associated with the "analytization" of the Lan- 
dau singularity; we shall revisit this issue in Sec. [7] and Appendix [D] Here the spectral 
density function pv([J?), as defined by Eq. ( 1B.6I) . stands for the timelike discontinuity of 
the coupling in the V-scheme, where 

A^ = A 2 e 5/3 , (3.19) 

where A 2 is defined in the MS scheme. As we shall see below (Eq. (I4.20p ) this is also the 
large-/3o limit of the "gluon bremsstrahlung" coupling. 



Dispersive representations of the Sudakov anomalous dimensions 

Eq. (13.181) with J- aing .(p 2 / m 2 , r) of (13 . 1 T[) constitute the dispersive representation of the 
Sudakov evolution kernel in the large-/3o limit; it involve two physical scales. We now wish 
to construct similar representations for the 'jet' and 'soft' Sudakov anomalous dimensions 
(13. 3p each involving just one scale. 

According to the leading-order result, Eq. (13. ip . ^ r sing .(/ i2 / m2 ? r ) must reach a finite 
limit for /i 2 — > 0, namely: 

-r fn ^ ln ( r ) , b i-d! 

^sin g .(0, r) = 1 — . (3.20) 

$As alternatives to the gluon mass regulator one can work in D-dimensions or use the Borel technique 
(see Appendix |B|) . 



14 



This, together with the functional form of T^Xl^ 2 / m 2 , r) summarized by (13.171) . implies 
that the two functions Tj and Ts are separately logarithmically divergent for /x 2 — > 0, 



Tj(e) 



ln(e) + h + 0(e) 



ln(e) + d x + O(e) . 



(3.21; 



Therefore, Eq. (13.181) . as written, cannot be split into the 'jet' and 'soft' contributions. 

An alternative representation of the kernel can be obtained using integration-by-parts; 
following [49] we define the timelike coupling ay ink (/z 2 ) through 



da^ nk (i2 2 ) 



< in V) = " 



dm? 



2 77T 



py(m 2 ), 



(3.22) 



obtaining 

R(m 2 , r) 



C 



R- 



large /3 



dfi 2 



Po Jo ^ 



C R 1 



df/ z 



Po r J /i 2 
where we use the usual convention, 



T 



3 



T s 



r 2 m 2 



(3.23) 



T(e,r) 



and similarly 



so 



Tj(e) = - 



Tj(e) 



^Tj(e); 
de 



dy <t ( \ 
—Tj{y); 

y 



dT(e, r) 
dlnl/e 



T s (e) 



T s (e) 



-4^(e) 
de 



(3.24) 

(3.25) 
(3.26) 



In the second line of Eq. (I3.23f) we inserted JF sing (/x 2 / 'm 2 , r) according to (13. 17ft . 

Infrared safety of ( 13.23[) implies that T sing ,(fJ> 2 / m 2 , r) should vanish for /i 2 — > 0. Obvi- 
ously this is not true for the separate terms Tj and Ts in the square brackets, which do 
not vanish for /z 2 — > but rather tend to a constant; according to (13.211) they obey 



Fj(0) = F s (0) = 1 



(3.27) 



Thus, the integral (13.231) cannot be simply split into two finite integrals one corresponding 
to the 'jet' contribution and one to the 'soft'. According to ( 13. 21) . however, such a split 
must be possible upon taking one logarithmic derivative. Indeed using ( I3.23P one gets 



1 dR(m 2 ,r) 



r dlnm 2 



CO 



and therefore 
dR(m, r) 



d\nm 2 



d^ 

large f3 P° JO ^ 
°° d/i 2 



p v (/j 2 )T Eing .(l^ 2 /'m 2 ,r) 



large a 



Po Jo I 1 ' 



-Pv{n 



T 



J 



rm* 



T s 



I'- 



r 2 m 2 



(3.28) 



(3.29) 
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which, in virtue of (13.271) . can be split into two finite integrals by subtracting = 1 and 

adding J\s(0) = 1 inside the square brackets. By comparing with (13. 2p . Eq. (13.291) implies 
that the 'jet' and 'soft' anomalous dimensions have the following dispersive representation 
in the large-/?o limit: 

,2 



1 r°° fin 2 r • ■ i 

^ 2 )L ge , =J Q J -j^Pv^ 2 ) [FAv 2 /k 2 )-Fj(0)\ , (3.30) 



and 



^ 2 )L g eA, =J o r^ Pv ^ 2) [^A 2 ) -JW)] , (3.31) 

where we used (I3.27p . Thus Tj and T s can be identified as the characteristic functions 
associated to the 'jet' and 'soft' Sudakov anomalous dimensions J and S, respectively. 
The corresponding representations in term of Oy ink (// 2 ) are: 

J(k%_, n =1 r^-a^(^j(^/k% (3.32) 

and 



llarge/3 ° A) Jo V 



^ k2 )L ge , =jJ o ^ar(M 2 )^ 2 A 2 ), (3-33) 



where 



Tj{e) = -e±Fj(e); T s {e) = -e± J^(e). (3.34) 

ae ae 

These representations of the Sudakov anomalous dimensions will be generalized in Sec. H] 
to full QCD. 

Dispersive representation of Sudakov anomalous dimensions as limits of the 
full physical evolution kernels 

Eq. (I3.30P and (13.311) can be obtained as r — > limits of the corresponding finite-r disper- 
sive representations of the physical evolution kernels in momentum space. We illustrate 
this statement (and identify the precise limit to be taken) with the examples of deep 
inelastic structure functions and the Drell-Yan cross section. This will also provide an 
illustration of the statement preceding Eq. (13. 18[) above. 

Deep inelastic structure functions 

Let us consider the large-/3o limit but finite x dispersive representation [49] of the physical 
evolution kernel. We note that at the level of the partonic calculation, the convolution 
on the r.h.s of (13.41) becomes trivial in the large-/9o limit since the 0(a s ) corrections to 
F 2 (z,Q 2 ) generate terms that are subleading by powers of 1 / /3 . Therefore, in this limit 
K dis (f 2 }(x,Q 2 ) is directly proportional to the partonic dF 2 (x, Q 2 )/d\nQ 2 , so 

(1 -x)ir DIs(F2) (x,Q 2 )| large/3o =-j?- J — py(/i 2 ) [(1 -x) (F(n 2 /Q 2 ,x) -F(0,x) 

(3.35) 
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where JF(/i 2 /Q 2 , x) is the standard notation for the characteristic function corresponding 
to F 2 (x, Q 2 ) (see Eq. (4.27) in [49]). In Eq. (13.351) we multiplied both sides by a factor of 
r = 1 — x in order to get a finite limit for the integrand for x — > 1 (see Eq. (13.451) below) , 
in agreement with the general result Eq. (13.1 7ft . Upon taking this limit under the integral 
with a fixed invariant mass W 2 = Q 2 (l — x)/x one obtains 



lim <^ (1 - x)K 



DIS(F; 



)(x,Q 2 



I large /3q 



x — > 1 
fixed W 2 



dfj 2 



Pv{^ 2 ) ^dis(^ 2 /W 2 ) - j- DIS (0) 



where we substituted 



^(/i 2 /^ 2 ) = I™ (1 - x)T^ 2 IQ 



(3.36) 



(3.37) 



x 



fixed W 2 

Eq. (I3.36P shows in particular that the specified limit of (1 — x)Kj 



DIS(F 2 ) 



>,Q 2 



I large /3q 



does 



exist. This (large-/?o) result is consistent with the more general result Eq. (13.91) (valid 
also at finite (3q). Moreover, Eq. (I3.9P allows to connect directly the momentum-space 
physical quantity K ms ^(x, Q 2 ) in the Sudakov limit with the physical Sudakov anomalous 
dimension, i.e. the exponentiation kernel. Comparing with Eq. (13.91) . we thus deduce 
Eq. (I3.30p . with the identification 

Fj(n 2 /k 2 )=F ms (n 2 /k 2 ). (3.38) 

The explicit result [57, 70] will be given in Eq. (13.451) below. 



Drell-Yan cross section 

In quite a similar way, we start from the large-/3o limit but finite r dispersive representation 
of the quark-antiquark partonic cross section g q q(r, Q 2 , fi 2 F ) in (13.101) . which was calculated 
in Ref. [49]. By taking a logarithmic derivative with respect to Q 2 we obtain the large-/?o 
result for the physical Drell-Yan evolution kernel (13.141) . 



[l-r)K DY (r,Q 2 )\ 



rKpv(^) (l-r)(j- DY (/. 2 /Q 2 ,r)-^ DY (0,r 



large/3 ° ft Jo v 



(3.39) 

where jF DY (/i 2 /Q 2 , r) is the characteristic function corresponding to g q g(r, Q 2 , fj, 2 ?) and 
where we multiplied by a factor of r = 1 — r. Next, by taking the r — > 1 limit under the 
integral while fixing the total energy radiated into the final state, E DY = Q{1 — r), we 
obtain: 



(l-r)K DY (r,Q 2 ) 



large j3 



C F f°° d[L 2 



~T / ^TPvW 2^ y (^ 2 /E 



T 1 

fixed E DY 



Po Jo V 



2-F DY (0) 



(3.40) 



17 



where we substituted 



2j- DY (/i 2 /^ Y ) = lim{(l - r^G^/QV)} 



(3.41) 



r -> 1 
fixed £/, 



and the limit can be shown [70] to exist, see Eq. (13.461) below. Comparing with Eq. (13.161) . 
we deduce Eq. (13.311) for S = S DY , with the identification: 



3.4 Characteristic functions of Sudakov anomalous dimensions: 



So far we have analyzed the properties of the physical Sudakov anomalous dimensions and 
their dispersive formulation in the large-/3o limit based on general considerations. Let us 
now summarize the results for the characteristic functions based on explicit calculations. 
These calculations have been usually done considering a specific distribution at finite r, 
computing the full characteristic function !F(e, r) by evaluating the squared matrix ele- 
ment with an off-shell gluon, and only at the end identifying the terms that are singular 
at r — > as demonstrated in Eqs. (13.371) and (13.411) above. There is however an alter- 
native: it was shown in [57] that there exists a systematic approximation — where one 
lightcone component of the gluon is taken small together with its transverse momentum 
squared and its virtuality — that allows a direct calculation of the r — > singular terms, 
namely Tj % s- The calculation of the soft characteristic function Ts can be simplified 
further, taking all the gluon momentum components to be small. This is the standard 
soft approximation albeit with an off-shell gluon. Here the hard partons that emit the 
radiation can be replaced by Wilson lines. The soft function can then be obtained through 
the renormalization of a corresponding Wilson-loop operator defined by the trajectories 
of the colored hard partons [32,54,85,86]. 

A compilation of the expressions for the characteristic functions in a few examples is 
given in Table [H The results for the derivative Tj^s-, defined by Eq. (I3.25P or (I3.26p . and 
the corresponding Borel function, defined by Eq. (16.70 below, are summarized in Table [2j 

Inclusive B decay and heavy-quark fragmentation 

Let us consider the characteristic function corresponding to the exponentiation kernel in in- 
clusive B decays. The simplest calculation is of the photon-energy spectrum c£T(m 2 , x)/dx 
in radiative B decays B — > X s j, where x = 2E 1 /rrib. This spectrum was computed in the 
large-/?o limit in [40]. The original calculation was done using the Borel representation 
of the dressed gluon propagator. In the dispersive formulation we obtain the following 



(3.42) 



results 
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process 



jet function 
(e.g. DIS) 



J 



B{y < I) 



i / n 3 1 1 2' 

Hv)- i + 2 y + l y 



process 



F s {z = e/r 2 ) 



B decay; HQ 
Fragmentation 



S x 



(JD 



ln(l + l/z) + 



1 + z 



Drell-Yan ; 
99 -» Higgs 



5 r; 



(z < 1/4) 



21ni±^H5-ln(z) 



e + e — >■ jets 
C parameter 
(r = c = (7/6) 



Sr. 



[z < 4) m(z) +9(z > 4) 



2 tanh 



1 (v / T 3 47i) - K*)' 



e + e — > jets 
Thrust 

( r = t = 1 - T) 



-6{z < l)In(z) 



Table 1: Results for the momentum-space characteristic functions of some inclusive dis- 
tributions. 



singular terms for r = 1 — x — > 0: 



•^"sing. (^5 ^) 



e < r 



— m(e/r) — 



le le 2 
2 r 4 r 2 



In (l + r 2 /e) 



1 + e/r 2 _ 



^7(e/r) 



^ QD (6/r2) 



(3.43) 



with e = fi 2 /m>l, where we have already split the result according to (13.1 7p and identified 
the two functions J r j(e/r) and jF5 QD (e/r 2 ). Equivalently, for the derivative we have: 



T (e r) = 

u sing. V, 5 / 



dlnl/e 



1 

r 



0(e < r) 



le 
2r 



2r 2 



+ 



e/r 2 



1 + e/r 2 (1 + e/r 2 ) 2 



^7 (e/r) 



^s QD (^ 2 ) 



(3.44) 



Here the superscript 'QD' stands for Quark Distribution, indicating that this particular 
soft function is related to the longitudinal momentum distribution of an off-shell b quark 
inside an on-shell b-quark state. The same result for the singular terms f!3.43j) was obtained 
for the triple differential rate in b — > X U W; here r is defined as the ratio between the two 
lightcone components p + /p~ of the X u jet, p ± — E ^ \p\, and the hard scale is p~ . 

Next we comment that both Tj{ejr) and ^(e/r 2 ) can be extracted from the calcu- 
lation of heavy-quark fragmentation in e + e~ annihilation, Eq. (54) in [30]. The fact that 
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process 




•^(y - e / r ) 




jet function 
(e.g. DIS) 


J 


%<1) [l-§-£] 


1/1 1 \ 
2\l-u l-u/2) 


process 






M s (u) 


B decay; HQ 
Fragmentation 




1 

1 (l + l^) 2 


71U 

1-u 

sin 7r-u 


Drell-Yan ; 
99 -> Higgs 




0(4z < 1) 

y/l-iz 


T 2 (l-n) 

r(i-2«) 


e + e~ — > jets 
C parameter 
(r = c = C/6) 


c 
<-> c 


1 > 4) 


r 2 (i + M ) 
r(i + 2u) 


e + e~ — > jets 
Thrust 

(r = t = 1 - T) 


<s 4 


1 - 6{z > 1) 


1 



Table 2: Results for the logarithmic derivatives of momentum-space characteristic func- 
tions of some inclusive distributions and the corresponding Borel functions (16. 7\i . 



^-5 (e/r 2 ) is the same for heavy quark decay and fragmentation follows from the general 
equality of the two Sudakov anomalous dimensions, proven in [32]. 

Finally, the quark distribution function at finite mass and arbitrary momentum frac- 
tion has been computed in the large-/3o limit in Ref. [32] using the Borel regularization. 
Repeating the calculation with a massive gluon and identifying the terms that are singular 
for r — > we recover ^ r ,s QD (e/r 2 ) of (I3.43p . 

The jet function 

As already mentioned in the introduction, the jet function appears in many different 
inclusive distributions where Sudakov logs emerge from a constraint on the invariant mass 
of a jet in the final state [57]. The simplest example is deep-inelastic structure functions, 
where J is the only source of Sudakov logarithms [5]. According to Eq. (13.371) . the 'jet' 
characteristic function is obtained upon considering the x — > 1 limit of the F 2 (x,Q 2 ) 
characteristic function (Eq. (4.27) in [49]) with a fixed invariant mass Q 2 (l — x)/x. The 
result is [57, 70]: 

J=>(e/(l-z)) = lim{(l -x)F(e,x)} 
x — > 1 



fixed e / ( 1 — x 
B(e<l-x 



(3.45) 



i 3 1 e 1 e ' 
- In - - + + -- 



x J 4 2 1 - x 4 (1 - x 
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where e = fi 2 /Q 2 . As anticipated, this result is identical to the jet-function contribution 
to various infrared safe observables including inclusive B decay (see [40]), Eq. (13.431) above 
with r = 1 — x, and a range of distributions in e + e~ annihilation into hadrons [57]: event- 
shape variables — e.g. Eq. (16) in [23] and Sec. 3 in [24], single particle inclusive cross 
section — Eq. (28) in [57], and heavy-quark production cross section. 

The Drell-Yan soft function 

The Drell-Yan (or Higgs) production case offers the basic example of a 'soft' function. 
Using in Eq. (13.411) above the explicit expression from Ref. [49] for the Drell-Yan charac- 
teristic function ^ r DY (e, r) with e = fi 2 /Q 2 , one obtains [70]: 



Ml 



r 



lim{(l-r)^ DY (e,r)} 



2 tanh 



fixed 



4e 



(3.46) 



6(4e/(l-r) 2 <l) 



(3.47) 



Upon defining z = e/(l — t) j 



T Sm {z) = 2 tanh~ 1 (y/l - Az) 6(z<l/A) = 



2 



Az 



-ln( 



z 



and therefore 



^ DY ( e /(l-T) 2 ) 



9(Ae/(l 



< I) 



46/(1 



(z < 1/4) 

(3.48) 

(3.49) 



t 



The same result has been obtained in Ref. [52] (see Eq. (4.6) there) and Ref. [57] (see Eq. 
(69) and the discussion following it) by considering directly the renormalon calculation 
in the soft approximation. Finally, note that in the Drell-Yan case the second derivative 
Fsby i s singular and therefore one must use (13.311) rather than (13.331) . 

Event -shape distributions 

Beyond the universal jet function discussed above, event-shape distributions in e + e~ an- 
nihilation are sensitive to large-angle soft gluon emission from the back-to-back recoiling 
quarks. The soft function strongly depends on the way the shape variable weighs the 
parton momenta. This is reflected in shape-variable dependent subleading Sudakov loga- 
rithms as well as power corrections. 

Before considering specific examples, a general comment is due concerning the not- 
completely-inclusive nature of event-shape variables: in contrast with the inclusive cross 
sections and decay spectra discussed above, event-shape variables do distinguish at some 
level between an off-shell gluon and the final-state particles to which it decays. Differences 
arise already in the large-/?o limit. The characteristic functions for event-shape variables 
are defined in the inclusive approximation, neglecting any such differences. 
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Let us first recall the result in the case of the thrust (or the heavy jet mass) which 
are particularly simple. Starting from the exponentiation kernel in the large-/3o limit (e.g. 
Eq. (20) in [23]) we have, in full analogy with Eq. ( 13.231) . 



C 



R(m 2 ,t) 



R- 



large fa 



Cr 



de 



2t 2t 2 



(e < t) 6(e > t 2 



(3.50) 



•^sing.^,*) 



where e = fi 2 /Q 2 , Q 2 is the center-of-mass energy and t — 1 — T, where T is the thrust 
variable. For thrust Cr = 2Cp. In order to extract the separate 'jet' and 'soft' contribu- 
tions to the characteristic function we use the identity: 

0(e < t) 0(e > t 2 ) = 0(e < t) - (l - 6(e > t 2 )) 

and note that the terms proportional to e/t or e 2 /t 2 in (13.501) contribute only at the e = t 
phase-space limit (they are power suppressed at the e = t 2 limit). This leads to the 
following decomposition of J r sing .(^,t) in (I3.50P 



•F sing. (^) t) 



6{e < t) 



1 



2t 2t 2 



(t>t z 



(3.51) 



V 



where we neglected regular contributions that are irrelevant for Sudakov resummation. 
Here we recognize the characteristic function of the jet from (I3.44p and (13.451) while for 
the soft contribution we have: 



^(e/t 2 ) = 1 - 9(e/t 2 > 1) = 9(e/t 2 < 1) 



(3.52) 



In a similar way we can write down the characteristic function of the C parameter 
using the results of Ref. [24]. Upon neglecting regular contributions, we obtain (see Eq. 
(3.5) in [24]): 



•£in g .(e,c) = 0(e < c) 9(e > Ac 2 



le_le 2 
2 c 2 c 2 



4c 2 /e 



l_l£_If! 
2 c 2 c 2 



VI - 4c 2 /e ( 1 + VI - 4c 2 /e 
(e > Ac 2 ) 



4c 2 /e 



^s c (e/c 2 ) 



(3.53) 



with c = C/6, where C is the conventional normalization of the C parameter (see [24]). 
Again we recognize the familiar jet function and find to be: 

* M = 1 - (3.54) 



1 - 4c 2 /e 

Note the similarity [54] with the case of the Drell-Yan soft function, Eq. (13.491) . Also here 
the second derivative jF Sc is singular, so one can only use (13 .3 1 j) . 
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4 From the large— /3q limit to the non-Abelian theory 



The dispersive technique is a convenient way to derive the all-order result in the large-/3o 
limit. But, in fact, it is much more than that: a dispersive representation, when it exists, 
summarizes the analytic structure of the observable in the complex momentum plane as 
well as its infrared sensitivity. Moreover, it is particularly well-suited to parametrizing 
power corrections by means of infrared-finite coupling. 

The purpose of the present section is to show that the dispersive representation of 
the Sudakov anomalous dimensions, Eqs. ( 13.30!) an d (13.311) . which emerges out of the 
calculation in the large-/5o limit, can be readily generalized to the full theory. In this 
generalization the characteristic functions, which have been computed analytically in the 
large-/3o limit, are kept fixed, while the V-scheme coupling ay iIlk , which was so far identified 
only to NLO in the large-/?o limit, is replaced by specific effective charges that are defined 
order-by-order in QCD and capture the non-Abelian nature of the interaction. After 
presenting the general formalism we will examine a few examples where corrections are 
known to the NNLO and beyond. The investigation of the evolution of these effective 
charges leads to interesting observations. 

4.1 A general dispersive representation of Sudakov anomalous 
dimensions 

The generalization of Eqs. (I3.30p and (13.311) (as well as (13.321) and (13.331) ) beyond the 
large-/5o limit is straightforward: 

"oo j, ,2 



Po Jo /i 2 L 

1 r ^VV)^ 2 A 2 ), 



(4.1) 



A) Jo /" 



and 



Po Jo L 



Po Jo P 2 
where, in full analogy with (I3.22p . 



(4.2) 



dahlia 2 ) _ f 00 dm 2 

pjAp 2 )= f; , ; *Ts(p 2 ) = - —pjA™ 2 ), (4-3) 

dmfi z ' J M 2 m z 

and, similarly to (IB. 61) . Pj,s(n 2 ) correspond to the timelike discontinuities of some "Eu- 
clidean" effective charges, originally defined for spacelike momenta: 



<"°o A, ,2 ,,2 lu2 /-oo 7 2 

'"' Pj.dir). (4.4) 



2 



/i 2 + k 2 
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where we ignored potential Landau singularities in a"(l; 2 ) that would violate the physical 
analytic properties of these effective charges (cf. Eq. (IB.5j) ). Of course, such singularities 
may appear at any given order in perturbation theory. We shall return to this issue in 
Sees. [O and [3 

Let us note in passing that while the Minkowskian representation of the Sudakov 
anomalous dimensions, Eqs. (14.11) and (I4.2p . can always be written, in general there are 
no equivalent Euclidean representations consisting of integrals over the effective charges 
a 1 j cl s (k 2 ). This fact reflects the timelike nature of these quantities. 

Obviously, we require that in the large-/?o limit ps(p 2 ) and pj(p 2 ) would both coincide 
with the discontinuity of the V-scheme coupling, pv(p 2 ) of Eq. ( IB. 61) . In Sec. 14.21 we will 
show, order-by-order in perturbation theory, that a generalization according to (14. II) and 
( 14. 2 p is indeed possible, and then further characterize these effective charges. We stress 
that the order-by-order analysis is entirely independent of the assumed analytic structure, 
as Landau singularities can only modify the dispersion relations (14. 4p by power-suppressed 
terms. 



4.2 From anomalous dimensions to effective charges: order— by- 
order relations 

We wish to show that the perturbative expansion of the anomalous dimensions, S and J ', 
and consequently their conventional decomposition in MS into the cusp A, B and V, 
Eq. (13. 3p . are effectively translated by (14. ip and (14.21) . order-by-order in perturbation 
theory, into an expansion of the corresponding effective charg es, a^(p 2 ). Thus, by 
virtue of the dispersion relations ( I4.4p . they uniquely translate into the expansion of the 
Euclidean effective charges, aFj cl s (p 2 ). 

Defining the expansion coefficients of aj^/i 2 ) in the MS scheme 



where a = /3o&f s (p 2 )/ 7! ', we have 



a* uc V) = 5 + J* a 2 + j*a 3 + jfa A + ■ 



(4.5) 



a™ ink (p 2 ) = a + j*a 2 + 



33 



a + sf a 2 + 



n 2 



S 3 ~ 



a 3 + 



a 3 + 



71 



5 
— < 
6 



a 4 + 



a 4 + 



(4.6) 



where we used (14. 4 p and the perturbative expansion of the running coupling in the MS 
scheme, where: 



a(zfi 2 ) = a(/i 2 ) - ln(z)a 2 (/i 2 ) + (ln 2 (^) -6\n(z)) a 3 (/i 2 ) 



+ -ln 3 (z) + -5\n 2 (z) -8 2 ln(z) a\p 2 ) + 



(4.7) 
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where the (3 function is given by 

da 



dm fi 2 



—a 



1 + 5 a + 5 2 a 2 + 



(4.8) 



with 5 = A/A, 2 and 5 n = f3 n /f3 n+1 . 

Next we note that Eqs. (14.11) and (14. 2p can be written as 

dy 



and 



Po Jo y 



dz 



A) 



(4.9) 



(4.10) 



Substituting p~7|) into (Ol) we get: 



^(y/i 2 ) = a(/i 2 ) + (jf - ln(y)) a V) 



+ 



7T 



j3 - - — + W( y ) -2HvW -sHv) 



a 3 (/x 2 ) + 



ln s (y)+ 3jf + ^ ln 2 (y) 



+ ( - 2jf (5 - 3^ + vr 2 - fc) ln(y) + jf - n 2 (jf + ^ 



a 4 (/i 2 ) + 



(4.11) 



and a similar expression for a^ ink (zfi 2 ) where j% are replaced by We observe that upon 
inserting these expansions into (14.91) and (I4.10p we would need to evaluate the log-moments 
of J-j{y) and j-s{z)i respectively, 



Jk 



Si 



o y 



00 dz 



(4.12) 
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The first few log-moments in some examples are given in Table [3] (see discussion below). 
Using Eq. (14.111) . Eq. (14. 9 j) yields the following expansions: 
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Po 
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+ {32 - Ji) a 2 + ( ji - y + J 2 - 2 Juf - J 1( J ) a 3 + 



[(a + a 2 a 2 + a 3 a 3 + ■ ■ • ) - (a 2 + 5a 3 + ■ ■ ■ ) (h + 2b 2 a + ••■)] ( 4 -!3) 



^l(a s (/i 2 ) + 



da(/i 2 ) dB(a s (ii 2 )) 
d In /i 2 da 



where we have identified the expansions of the anomalous dimensions A(a s (fi 2 )) and 
B{a s {fi 2 )) of Eqs. (J23D and fl23D) , respectively. At order C(a 2 ) we find 6i = Ji- Eq. (14331 
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Table 3: Summary of results for the first few log-moments of the characteristic functions, 
Jfc and Sk with k = to 5, of some inclusive distributions (cf. Table [2]). The log-moments 
are defined in Eq. (14. 12p . They can be most easily extracted using (16. 7p by expanding the 
Borel functions in Table [2j 



yields the following expressions for the coefficients for the Euclidean effective charge a?j ' 
in gU): 



32 = a 2 



jf = a 3 - J 2 + 2b x a 2 + -vr 2 - 2b 2 



(4.14) 



31 = a 4 + 66 2 a 2 + (2a 2 <5 - 3 J 2 - Qb 2 + 3a 3 ) b x + ^-3a 2 - -5 J J 2 + J 3 
5 

+ -ir 2 5 - 3b 3 + vr 2 a 2 - 25b 2 , 
6 

where the log-moments Jk are defined in (14.121) and a n and b n are the MS coefficients of 
(E5D and (EHUD, respectively. 
Similarly ( 14.101) yields: 



/3q 



7l~ 



+ (sf - Si) a 2 +(s%- — + S 2 - 2S lS % -S 1 6)a 3 + - 



C R 



A(a s (n 2 ) 



[(a + a 2 a 2 + a 3 a 3 H ) - (a 2 + 5a 3 H ) (di + 2d 2 a H )] (4.15) 

dd(/j, 2 ) dV(a s (/j, 2 )) 



cHn/i 2 



where the expansions of the anomalous dimensions A(a s (/i 2 )) and V(a s (fj J 2 )) of Eqs. (I2.9P 
and ( 12. lip , respectively, have been identified. At order 0(d 2 ) we find d\ = S±. Eq. ( 14. 15ft 
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yields the following coefficients for af d in g3D 

E 

s 2 = a 2 

sf = a 3 - S 2 + 2d x a 2 + \v; 2 - 2d 2 

o 



5 



(4.16) 



s 4 = a 4 + 6rf 1 a 2 + (2a 2 5 - 35*2 - 6d 2 + 3a 3 ) c?i + ^-3a 2 - ~8J S 2 + 5 3 
5 

+ -n 2 5 - 3rf 3 + n 2 a 2 - 25d 2 , 
6 

where the process-dependent log-moments Sk are defined in (I4.12p and a n and the process- 
dependent d n are the MS coefficients of (12.91) and (12. lip , respectively. 

Obviously, such an identification can be done to arbitrarily high order: according to 
(14. ip and (14.21) . the information contained in J and S at a given order in the perturbative 
expansion is sufficient to determine a^j cl and a| ucl , as well as aj mk and a^ ink , to that order. 

Finally, note that following the identification we made above we can write the small-e 
expansion of the characteristic functions (I3.2ip as 

Tj{e) = - ln(e) + h + 0(e) = - ln(e) + f° ^ ln(y) tj(y) + 0(e) (4.17) 

jo y 

F s (e) = - He) + S 1 + 0(e) = - ln(e) + / — ln(^) F s (z) + 0(e) (4.18) 

Jo z 

where the constant terms b\ = J\ and d\ = Si were fixed by the first log-moment of the 
characteristic function itself. 

Cusps and NLO universality of the "gluon bremsstrahlung" coupling 

We observe that the effective charges Oj 1 ^^ 2 ) are closely related to the cusp anomalous 
dimension [65,66]: 

A ^ <k or ^ ?* aT ^[d + a 2 d 2 + ■■•]; (4.19) 

Po Po Po 

their expansions start deviating [68,69] only at the NNLO. We emphasize that the NLO 
universality of the effective charges holds despite the fact that the source of the soft radia- 
tion can be quite different: it applies to the entire class of inclusive distributions discussed 
in this paper, including B decay spectra, event-shape distributions, deep inelastic struc- 
ture functions, single-particle inclusive cross sections, Drell-Yan production and heavy 
quark fragmentation. To NLL accuracy, the only way by which the non-Abelian nature 
of the interaction affects the exponent is through the cusp anomalous dimension (see 
e.g. [32,85,86]), and therefore this effect is independent of the details of the configuration 
of color sources emitting the radiation. 

This NLO universality explains and justifies a posteriori some results obtained in the 
past in several occasions [12,22,49]. In particular, in [12] it was proposed to absorb 
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the NLO term into the definition of the coupling, to define the "gluon bremsstrahlung" 
coupling. This can be formulated as a rescaling of the dimensional transmutation scale A 
as follows (see also Eq. (60) in [12]): 

A L = A 2 exp {af } =A 2 exp{j; + ^Q-£)}, (4.20) 

where the subscript GB stands for "gluon bremsstrahlung" and where A 2 is defined in MS. 

Thanks to this NLO universality, in the original DGE paper [22] it was possible to 
generalize the exponentiation of a single dressed gluon to full next-to-leading logarith- 
mic accuracy, without having to introduce two separate effective charges. At this level, 
differences between observables, which affects subleading logarithms as well as power cor- 
rections, are encoded only in the different functional form of the characteristic function. 
Beyond the NLO the 'soft' and the 'jet' effective charges start differing from the cusp 
anomalous dimension, which also becomes scheme dependent. 

Note that differences with respect to the cusp anomalous dimension appear already in 
the large-/?o limit. Nevertheless, we expect that similarly to the cusp anomalous dimen- 
sion, the effective charges a^j cl 0.') and would be reuormalon freel Thus, the 
expansions (14.51) should not develop large subleading corrections. Absence of renormalons 
in cijr ucl (/i 2 ) and a| ucl (/x 2 ) implies that also j7"(/i 2 ) and S(fi 2 ) are renormalon- free, since the 
small fi 2 /k 2 expansion of the characteristic functions in (14.11) and (14.21) does not give rise 
to any non-analytic terms. These relations will be further discussed in Sec. 16.11 using the 
Borel formulation. 



Examples 

Let us now use the explicit expressions for the characteristic functions in Sec. 13.41 to- 
gether with the known perturbative expansions of the Sudakov anomalous dimensions, 
Eqs. (I2.9p . (I2.10p and (I2.1ip with the coefficients in Appendix [S] to derive the correspond- 
ing expansions of the effective charges a|j cl and a| ucl . 

The jet function 

Starting with the jet function case, from (I3.44p we get: 
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Tj{y) = 6(y<l)y i - - u 



(4.21) 



so using (14.91) we obtain 

1 rdy 
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o Jo 
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aT k (yv 2 )FAy) 
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dy a- nk (2//i 2 ) y + 



(4.22) 



§It should be emphasized that absence of renormalons in aj- ucl (/x 2 ) and a| ucl (/x 2 ) is obvious in the 
large-/?o limit, where these functions coincide with the one-loop coupling (|4.39p . but it is far from obvious 
beyond this limit. 
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The corresponding log-moments defined in (14.121) are: 



Jk 



^-hx k {y)Tj{y) 
o V 



dy \n k (y) [y+- 



-l) k k\ (1 + 2~ k )/2 . (4.23) 



Note that these coefficients increase factorially. The first few numerical values are sum- 
marized in Table [3j 

Using now (14. 14[) with the MS b n in Appendix [A] we get: 
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(4.24) 

where the coefficients a n , which are known in full to three-loop order (see Appendix |A|) 
have not been substituted. As with the anomalous dimension J itself, the only missing 
ingredient in ajj cl (/i 2 ) at four-loop order (N 3 LO) is the four-loop coefficient of the cusp 
anomalous dimension, a^. the combination CrcFj c1 / (3q — A is known to this order. 

The fact that the four-loop coefficient is not known in full is rather unimportant in 
practice: the series for the cusp anomalous dimension A in MS converges very well, and 
a 4 can be reliably estimated using Pade approximants [87] as discussed in Sec. 4 in [7]. 
For example, for N c = 3 and Nf = 4 we get: 



A 



A) 



a + a 2 ■ 0.962 + a 3 • 1.159 + a 4 



1.341 [1,2] Pade 
1.397 [2, 1] Pade 



+ 



(4.25) 
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where the difference between the two Pade approximant predictions can be used as a rough 
measure of uncertainty. 

Using these coefficients of A in (14.241) we getQ: 

a T = a + a* • 0.962 - a* • 1.019 - * ■ { V?' ^ f' *j " » } + • . • . 

J [3.165 using [2. 1J Pade for cu m (I4.25P J 

(4.26) 

Note that this expansion converges very well, better than the corresponding Minkowskian 
coupling: 



a™ ink = a + a 2 ■ 0.962 - a 3 ■ 4.309 - a 4 



18.798 using [1,2] Pade for a 4 in (ji^SD 
18.742 using [2, 1] Pade for a 4 in (T4T251) 



(4.27) 



and the physical anomalous dimension J itself: 



,<?o J = a + a 2 ■ 1.712 - # . 1.061 - a* . { £™ ^ £• *] ta } + . . . , 

[ 17.549 using [2, 1J Pade tor a 4 m ( 14.251) J 

(4.28) 

which both contain large 7r 2 terms (see (14.61) ) owing to the analytic continuation to the 
timelike region. It should be emphasized that (at least in the large-/?o limit) none of the 
three has renormalons, despite the fact that the log-moments increase factorially. The 
most convenient way to see this is by considering the relations between these quantities in 
terms of their Borel representations, see Sec. 16.11 below. 

The Quark Distribution function (B decay) 

Consider next the soft function associated with the Quark Distribution in an on-shell 
heavy quark. From (I3.44[) we get: 

^ ) = 1 -(TTW ^ ^ (2)= (TTIF (429) 

Using (14.101) we get 

Sqo ^ 2)= WoL t^^^W = jj ^CWji^- ( 4 - 3 °) 

The corresponding log-moments Sk are given by 



ST=I ^ln k (z)F SQO (z)= I dz\n k (z)-—-. (4.31) 



Their first few values are summarised in Table [31 

^Note in contrast with (|4. 25[) . a direct Pade approximant prediction of the N 3 LO coefficient in aj ucl is 
unreliable: the [1,2] and [2,1] approximants differ by much. A direct prediction for (14.27)1 or (14.28[) below 
is even worse owing to the large ir 2 terms related to the analytic continuation. 
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Using (14.141) we get the following results for coefficients in Eq. (14.51) : 
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Numerically, we find again very good convergence in the first few orders; for N c = 3 and 
Nf = 4 we get: 



<5qd 



a + a 2 ■ 0.962 - a 3 • 0.799 + 



(4.33) 



The Drell-Yan soft function 

Consider now the soft function associated with Drell-Yan production or Higgs production 
through gluon-gluon fusion. The characteristic function is given by (I3.49p . Here the second 
derivative is singular, but using (I3.3ip one can verify that the relations in Eq. (14.161) hold. 
Then, using (16.71) to extract the log-moments Sk (see Table |3]) one readily obtains: 

°|£(m 2 ) = s + s 2^ 2 + s 3^ 3 + s i « 4 + • • • 



28 C A /8 7 . 
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240 
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7 49 \ 3 / 11 77 
18 + 32 C3 J CA + ri8 + 32 



+ ^C 3 ) C F C A 2 



a 4 + 



(4.34) 



where we used the known coefficients ( 1A.3I) . Using (1A.1I) . Eq. (14.341) yields the following 
expansion for N c = 3 and Nf = 4: 



o§£ = a + a 2 ■ 0.962 - a A • 0.508 - a 



^3 



9.07 using[l,2]Padefora 4 in(ji25D 
9.01 using [2, 1] Pade for a 4 in fl4T25l) 



+ 



(4.35) 
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4.3 Renormalizat ion— group evolution of the Sudakov effective 
charges 

In order to evaluate the anomalous dimensions entering the Sudakov exponent one needs 
to perform integrals over the discontinuity of the running coupling in (14.11) and (14.21) . The 
scale dependence of the relevant Euclidean couplings can be computed directly [63] by 
solving the appropriate renormalization-group equations. Upon taking the derivative of 
the effective charges a^ ucl (/x 2 ) and a| ucl (/x 2 ) with respect to ln/i 2 and expressing the result 
in terms of the effective charge itself one arrives at the following renormalization-group 
equations: 



din n 2 



Si = Sr + 3i-3%S-(jZ) 2 

S( = 5T S + 2 j? + 4 (if ) 3 + 5 (jf) 2 - 6jfjf - 2j*%* 



(4.36) 



= - (*v» 



1 + arWS + (aTW) 2 S S 2 + (a?V)) S *? + 
SS = Sr + si-s*S-(s*) 2 

Si = <5™ + 2af + 4 3 + 5 ( S f) 2 - 6*f S f - 2siSf s 



(4.37) 



These equations can be truncated at any given order, and integrated. The first ap- 
proximation is obtained setting 5 = 5 n = which is also the result of the large-/3o limit, 

«rv)L g « ft = wiu, = rrjfc • (438) 

Upon using the universal NLO correction of (14.191) to fix the initial condition one obtains: 
af c \fi 2 )\ . = a E 7 uc V)l , =- J——, (4.39) 

15 V.^ / 1 one-loop J \t~ / \ one-loop J n ^2 j J^2 ' v ' 

where A GB is defined in Eq. (14.20}) . At the next truncation order (5^0 but £2,3... = 0) the 
effective charges aj cl (/z 2 ) and a| ucl (/x 2 ) still coincide. The analytic solution can be written 
in terms of the Lambert W function [72,73], facilitating exact analytic continuation to the 
time-like axis. 

Beyond this order aj cl (/i 2 ) and a| ucl (/x 2 ) start differing from each other. It is only at 
this level that the process-dependent nature of soft gluons radiation reveals itself. Indeed 
in the large-/?o limit 5 n = 0. The evolution equations (I4.36P and (I4.37P directly reflect the 
non-Abelian nature of the interaction. It therefore becomes particularly important to 
study the higher-order corrections to these equations. 

Here comes a very interesting observation: in contrast with the MS coupling and with 
other physical Euclidean effective charges [72], the higher-order coefficients &f 3 are not 



32 



all positive. Most importantly 8% is negative in all the known examples (see below). 
This opens up a possibility that an infrared fixed-point would be realized already at the 
perturbative level. In this case Landau singularities may not appear in the entire first sheet 
of the complex momentum plane, leading to a causal analyticity structure [72,74], where 
the entire complex momentum plane is mapped into a compact region in the complex 
coupling plane. A Landau singularity does appear in the one-loop coupling (14.391) . in 
the large-/3o limit, as well as in the two-loop coupling with any realistic number of light 
flavors. Its absence would therefore be a direct consequence of the non-Abelian nature of 
soft-gluon interaction, which may be different in different processes. 

If this scenario is realized the perturbative solution for the effective charge can genuinely 
be a good approximation to the full, non-perturbative effective charge down to the infrared 
limit. In this case the Sudakov factor computed as a dispersive integral will be in very 
good control, far beyond what can be achieved in a fixed-logarithmic-accuracy approach. 

Of course, in reality we know, at best, just 82 and 03 and therefore it is hard to 
make any firm conclusions about the size and perturbative stability of the fixed point. 
Obviously, a fixed point can always be washed out by sufficiently large subleading correc- 
tions. Nevertheless, some of the examples considered below (in particular the Drell-Yan 
case) are quite suggestive of this scenario. 

It is important to stress that the infrared finite coupling scenario discussed here has 
nothing to do with Analytic Perturbation Theory (APT) coupling [88], which becomes 
finite only owing to the "analytization" of the Landau singularities by imposing the disper- 
sion relation [33,38,88]. This model implicitly assumes the existence of power corrections 
that are not of infrared origin. See further discussion of this issue in Sec. [7] and Appendix 

El 

To examine the effective-charge beta function in Eqs. (14.361) and (14.371) it is useful 
to introduce a decomposition of the coefficients into powers of (3q (eliminating the Nf 
dependence in favor of po = j^Ca — |^/)- To this end we define first the decomposition 
of the coefficients of the anomalous dimensions bi and di in the MS scheme, Eqs. ( 12. 91) . 
(I2TTU]) and (|2TTTp . respectively: 

j i— 1 ^ 4—1 ^ 4—1 

a i = -7^i^2 a i,k^ k ; k = —i^h^- d i = — [ZT J2 d i,kPo k (4-40) 
and the MS beta function 

i 

j=0 

with (5 4 MS = p\ MS /pV +1 - In what follows we omit the superscript MS, and simply use 
when referring to this scheme. Using similar notation we also define the decomposition of 
the coefficients of the effective-charge beta function in Eqs. (I4.36P and (14.371) : 

? = |i = iE € /V; € <V , (4.42) 
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where the fact that the sum over j goes only up to j = i, as in the the MS scheme (14.411) . 
rather than j — i + 1 (as occurs in most physical effective charges) is a reflection of the 
fact that in the large-/3o limit the Sudakov effective charges are simply related to the MS 
coupling, see (14.381) above. 

Considering the first scheme-dependent coefficient, 

rj ^2 X J , $2,2 $2,1 $2,0 X S P2 X S , ^2,2 $2,1 ^2,0 , Q v 

$2 =-73= $2,3 +~jT + + ^3 i $2=^3 = $2,3 +~jT + ^2 + ^3 ( 4 ' 43 ) 

Po Po Po Po Po Po Po Po 

=0 =0 

we observe, as anticipated, that 5 2 ^ s = 0. We also find that S 20 = $20 an d $21 = $21 are 
universal, i.e. they depend only on the coefficients of the cusp anomalous dimension and 
the MS beta function (they do not depend on the details of the soft function considered) 

7 S 

while S 22 are process dependent, depending on bi = J\ and 62 or on d\ = S\ and g?2, 
respectively. The expressions are: 

^2,0 = ~~ a 2,0p\,0 + P2,0 

$2,1 = ~ a 2,0 + P2,l ~ a 2,lp\,0 — a 2,0p\,l + a 3,0 

$2,2 = -2a 2 ,oa2,i + 2S , ia 2 ,o + «3,i - o 2 ,iA,i + £2,2 - 2d 2 , ( ' 4 ' 44 ' ) 

7T 2 

$2,3 = -^2 + 2S , ia 2 ,i + — - 2d 2 ,i - a\ x + a 3i2 = , 
and similarly for Sfj, with the obvious replacement of the coefficients dk by the corre- 

7 S 

sponding 6 fc and S*. by J*.. Note that the vanishing of 5 2 ^ can be verified by taking the 
large-/3o limit of 5 2 ' in Eqs. ( 14.361) and ( 14. 3 7ft and using jf and sf of Eqs. (14.14p and 
(14.16p . respectively. 
At the next order, 

sJ xJ xJ xJ XS xS 

xJ _ xJ ,^3 ^2 ^1 ^0. ,5 _ ,5 ,^3 ^ Al «3j0 f4 4 , 

v - v ^' Po Po PO PO ^-v-' Po Po Po Po 

=0 =0 

so we find, as expected, that 5f^ s = 0. We also observe that &f = <5f is universal, 

/ 7 S 7 S 7 S\ 

while the other coefficients [o^[ and £32 and £33 ) are process dependent, depending 
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on d\ = Si, d 2 and <i 3 as well as on S 2 . The general expressions are: 

^3,0 = — 2/32,0^2,0 + /3l,0°2,0 + 03,0 

<jf tl = 2a 4 , + ^i.ia^o + 4a2, + /? 3 ,i - 2/92,002,1 

+ 2/3i,oa 2 ,oci2,i + 4a 2 ,o/3i,o<Si - 4/3i i0 d 2 ,o - 2/3 2 ,ia2, - 6a 3 , a 2 ,o 

5 f,2 = ( 6a 3,o - 12a| + 4a 2 , /3i,i + 4a 2 ,i/?i,o) Si + 2a 4 ,i 

571-2 

+ /?3,2 - 6a 3 , a2,i + 12a2, a 2 ,i + -3-^1,0 - 6d 3 ,o 

- 2/3 2 , 1^2,1 + /3l,0 a 2 1 — 2/^2,2^2,0 

(4.46) 

- 6a 2 , a 3 ,i - 4/3i, d 2 ,i - 4/3 1 ,id 2 ,o + 2/3i, 1 a 2 , a 2j i + 12a 2 ,od 2 ,o - 5/3 h0 S 2 
5f 3 = 12S , ia 2 ,o + (-24a 2 , a2,i + 6a 3 ,i - 12d 2 ,o + 4a 2 ,i/3i,i) Si 

+ /3 3 , 3 + 12a 2)0 a 2)1 ~ 2/92,2^2,1 - 6d 3 ,i + /3i,ia 2)1 - 5Pi,iS 2 - 4/3i,id 2 ,i 

571-2 

+ 2a 4 , 2 + Y2a 2 fld 2 ,i - 6a 2 ,oa 3 ,2 + — ^— /?i,i + 12a 2 ,id 2 ,o - 6a 2jl a 3 ,i 

5f j4 = 125^2,1 + (-12a^ + 6a 3 , 2 - 6S 2 - 12d 2 ,i) Si + 25 3 

- 6a 2 ,ia 3 , 2 - 6d 3 ,2 + 2a 4 , 3 + 4a 21 + 12a 2 ,ic/ 2 ,i = 

and similarly for 8fj with the obvious replacements. Again, the vanishing of Sf^ s can be 
verified by taking the large-/?o limit of 5% ,s in Eqs. (I4.36P and (14.371) and using jf and sf 
in Eqs. ( 14. 14ft and ( 14. 16[) . respectively. 

We comment that the above universality of 52.0; ^2,1 and 5 3 ,o is directly related to the 
universality of the Banks-Zaks fixed point [68,89], which will be proven (to all orders) in 
Sec. 4.4 below. This relation follows immediately from Eqs. (3.3) and (3.6) of [67], taking 
into account the universality of the Banks-Zaks critical exponent. 



Three— loop examples 

Let us now examine the coefficients of the effective-charge f3 functions for the examples 
considered in Sec. 14.21 Using Eqs. (14.441) and (I4.44|) with the explicit coefficients in Ap- 
pendix |A] and Table [3] (or, alternatively, using in (14.361) and (I4.37P the explicit coefficients 
fl4T24j) . (|4T32|) and (IQijl ) together with the MS /3-function coefficients, we obtain the fol- 
lowing results for the three-loop coefficients for the effective-charge j3 functions: for the 
jet function (I4.24p we find 



fi J l) J ti J 

°2,2 °2,1 °2,0 
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with 



^ = (s - 2 6 + r J ^ + ( is ' r) Cf 

sil = + + z? _ ha c/ + /as + h + vj CfCa _ £ c/ 

2jl V 120 48 96 4 7 V 96 4 s3 16 J 32 F 

2 '° V 512 192 / V 64 192 J 128 

(4.48) 

Another constraint on the jet Sudakov anomalous dimension is the color structure. In 
particular, for Nf = we have: 

o7\ f 121 121 2 121 A^2^ A? 29 121 . 21 2 11 4 \ 3 

#fL-n = - 7^ + "TTtCs C A 2 C F +[—- — ( 3 + 77^ + 77777^ ^ , 



^=0-^768-576* + ^ C 7 A F+ Vll52 - ^ 2 " l^ 77 + I 449 ' 

(4.49) 

where the Cp 3 and the C A Cp 2 components vanishes. In (14.471) the latter implies the 
following relation between the coefficients of (14.481) : 
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$?f F ) + <f F ° A) = • ( 4 -50) 

This relation involves only the universal components of &f that coincide with the corre- 
sponding ones in 8% , so it holds for any Sudakov anomalous dimension. 

Let us turn now to the soft anomalous dimensions. For the quark-distribution function 
( EQ2J) we find 

j-Sq D cSqo rS QD 

Po Po Po 

with 

'23 \ „ /97 



SiT 3Cs J Cf + 1,32 " Cs J Ca (452) 

and S^f = and 82 q D = &f , given in Eq. (I4.48P above. For the Drell-Yan soft function 
(Q4j) we find' 

xSdy A^ dy xSdy 

tf" = T + 3T + %r. (453) 

Po Po Po 

with 

^ Y = (f " 3Cs) C F + gcU (4.54) 
and 8^ Y = and 5*> Y = <Jf , given in Eq. (14.481) above. 
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The constraint on the color structure of the soft anomalous dimensions is more stringent 
than for the jet: they are maximally non-Abelian. In particular, for Nf = we have: 

(4.55) 

H2 Nf =o \32 288 36 1440 ) 

where the Cp 3 , the Cp Ca and the CfCa 2 components all vanish. Thus, beyond (14.501) . 
which is obviously realized, there is an additional relation, namely the Cp component of 
#f 2 is universal: 

S (c F ) _ 144 j, S {c F cl) 11,j,s(c a c f ) _ 23 

' 2 _ _ 121 2,0 ~ 11 2,1 _ ~8~ (4.i>t>j 

This general property is of course realized in the examples above, Eqs. (14.521) and (14.541) . 
Thus, at this order, only the Ca component of <5f 2 distinguishes between different soft 
functions. 

Incidentally, the numerical values of &f an d <5f are also not far; for example for 
N c = 3 and N f = 4 (where 5 = 0.7392) one gets: <Jf ~ -1.954, £f QD ~ -1.734, and 
5^ DY ~ -1.443. 

As anticipated, the numerical values of 5% and 5% are all negative, allowing for a 
perturbative infrared fixed point. It thus becomes interesting to examine higher-order 
corrections to the effective-charge (3 functions. 



5dy 



JV/=0 



523 17 , 

1 7T 2 

288 288 
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Four— loop examples 

In case of the jet function and the soft function associated with Drell-Yan production we 
do have sufficient information to determine the corresponding four-loop coefficients 63 of 
Eq. (14.461) using the Pade approximant predictions for 04 in (14.251) . For N c = 3 and Nf = 4 
we get the following effective-charge (3 functions: 

day i {fi 2 ) 



dlnfi 2 



(arv)r x 



1 + 0.739 • c4 ucl (/i 2 ) 



(4.57) 



and 

<fag£(/* 2 ) 
dlia.fi 2 



2\\ 2 



X 



1 + 0.739 . a |£ (fi 2 ) 



(4.58) 



"The corresponding value quoted in [68] is incorrect. 
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In case of the jet function effective charge (14.571) the positive £3 coefficient does not support 
the fixed-point scenario; a more definite conclusion requires higher orders. However, in 
the case of the soft Drell-Yan effective charge #3 is negative and sufficiently large to bring 
the fixed-point value to the perturbative regime: at this truncation order we get a zero for 
the j3 function at as DY ~ 0.46 corresponding to af DY ~ 0.7. Of course, since there is no 
perturbative stability at this order, higher-order corrections may well change this infrared 
fixed-point value. Further insight on this issue will be gained in the next section. 

4.4 The universal infrared limit of Sudakov effective charges 

If we assume that the Sudakov effective charges do indeed admit a finite infrared limit, it 
becomes natural to ask what this limit is, and further, to what extent it is universal. The 
most natural tool to address these questions is the Banks-Zaks expansion [72-74,90-93]. 

Banks-Zaks fixed point 

Let us briefly recall some terminology. A Banks-Zaks fixed point [90] is a non-trivial 
infrared fixed point, (3 (a(/z 2 — > 0)) = 0, occurring in asymptotically free gauge theories 
with a sufficient amount of matter. In QCD, such a conformal infrared limit characterizes 
the theory with Nf/N c 33/2 where the beta function 

p (cV)) = ^[^J 71 = -A) {a s ^ 2 )/n) 2 ~fa {a s ^ 2 )/nf-fa (a s (/i 2 )/vr) 4 + - • • (4.59) 

has fa > and fa < 0, so it is negative for a vanishingly small a s (fi 2 ) admitting asymp- 
totic freedom, but then changes sign again at a s (fi 2 ) ~ — fa/ fa, implying that theory is 
conformal in the infrared limit, where the coupling saturates at a s (fi 2 — > 0) ~ —fa/ fa. 
In the fa —>■ limit the infrared physics is fully under control of perturbation theory. In 
real-world QCD both fa and fa are positive, and for most physical effective charges so 
are the higher-order coefficients fa, fa etc.; then there is no infrared fixed point. 

The Banks-Zaks expansion [72-74, 90-93] amounts to expressing the infrared fixed- 
point value of the coupling as a systematic expansion in (positive) powers of fa. To 
construct the expansion one eliminates the Nf dependence in favor of fa = %Ca — \Nf, 
and then solves the equation (3 (a) = (the tilde denotes a particular renormalization 
scheme) order-by-order in fa, using the fact that in the formal fa — > limit the infrared 
coupling itself is of the order of fa. Using the notation of the previous section for the 
decomposition of the effective-charge beta function coefficients, the resulting expansion 
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takes the form (see e.g. Sec. 3.1 in [72]): 



a s (/i 2 -> 0) 



7T 



1. 1 



(ho 



2 ft,o + 3A,ift,o | 24% 

A,0 /?% 



+ 



Aj,2 /9i,o + /9?i + 3/3 a>1 A, 1 + ^3,1 - 



4 gi, i P 3 , o + 6 g i o + 4 gg, fo. 



+ 



10/3 



1. 1 



+ 5/3 2 , /33,( 



5 A 3 



2.0 



+ 



(4.60) 



where b BZ = — Po/Pi,o- This expansion converges well for an appropriate number of light 
quark flavors (Nf/N c ^ 33/2) where a fixed point exists (and is realized perturbatively) 
while for a realistic number of light flavors, its convergence strongly depends on the ob- 
servable considered [72]; this is consistent with the expectation that real-world QCD does 
not have a conformal infrared limit. Physical quantities may have a finite infrared limit 
also in the confining phase, which is not driven by a conformal fixed point and is usu- 
ally inaccessible to perturbation theory. Yet, it is possible that certain physical effective 
charges can be described by perturbation theory down to the infrared limit and then the 
Banks-Zaks expansion provides a natural tool to compute it. 



The universal infrared limit of Sudakov effective charges 

In general the infrared limit of the coupling is scheme dependent. This is true also for 
physical effective charges and it is reflected in observable-dependent coefficients of the 
Banks-Zaks expansion starting at NLO [72,91]. Quite remarkably, and in sharp contrast 
with the general situation, Sudakov effective charges have a common infrared limit that is 
uniquely determined by the cusp anomalous dimension; as we shall see, their Banks-Zaks 
expansion is universal to all orders in perturbation theory. 

This type of universality has been first observed in Ref. [89] where it was noted that 
the Banks-Zaks expansion of the Sudakov effective charges relevant to deep inelastic struc- 
ture functions and Drell-Yan production coincide up to the N 3 LO. The issue was further 
developed in Ref. [68] where a connection was made with the infrared limit of the cusp 
anomalous dimension. Here we give a general prooQ of the above statements. 

To prove the universality of the infrared limit note first that the relations between 
the physical Sudakov anomalous dimensions and the MS ones (13.31) involve a momentum 
derivative, thus the f3 function itself. Upon assuming that all effective charges involved 
admit an infrared fixed point, Eq. (I3.3P immediately implies that in this limit the term 
containing the f3 function vanishes, so [68]: 

C R J(fi 2 -»• 0) = C R S(fi 2 -> 0) = A(fi 2 -> 0). (4.61) 

**A different argument, which also applies to the more general situation where the infrared fixed points 
are of a non-perturbative origin and not necessarily related to the cusp anomalous dimension (which may 
have no fixed point at the non-perturbative level) is given in [68]: see Eq. (27) there, where instead a 
relation to the quark form factor is established. 
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The next observation is that the dispersion relations (14.11) and (14.21) necessarily relate the 
infrared limit of the Sudakov anomalous dimension with that of the corresponding Sudakov 
effective charge: 



3^ - 0) = i - 0) = 1 a^V -> 0); 

Po Po 

5(/, 2 - 0) = 1 ar k (/i 2 - 0) = i- aT (M 2 - 0) , 

Po Po 



(4.62) 



where the equality of the Minkowskian and Euclidean effective charges [72] follows directly 
form (I4.4p . It is worthwhile emphasizing that this equality holds despite the fact that the 
function of the Minkowskian effective charge differs from the corresponding Euclidean one 
by large 7r 2 terms at any order (at three-loops and beyond). The Banks-Zaks expansion 
of a Minkowskian effective charge is identical to the corresponding Euclidean one (Sec. 
3.3 in [72]). 

Using in (I4.60P the results for the effective-charge beta function of the Sudakov effective 
charges, Eqs. fT4~44l and (IPfi]) . we find, as implied by fT4~611 and (EE2), that all the 
process-dependent coefficients drop out: at each order in the Banks-Zaks expansion only 
the cusp anomalous dimension and the QCD beta function coefficients appear. The result 
for the first few orders is: 
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1,0 



01, ( 



(4.63) 



where the coefficients of the cusp anomalous dimension Ojj and the QCD beta function 
0jj are defined in (I4.4(jp and (I4.4ip . respectively. Both these quantities are defined in the 
MS scheme, but this scheme dependence cancels at each order in the expansion, as both 
the l.h.s and the expansion parameter are physical quantities. 
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Finally, for N c = 3 QCD we find the following numerical values: 
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4750031643031 75857683 



4955688493 5335 



422055936 



^56 



C3 J " 



1083839643648 3297312 



0.3308 b 2 BZ + 6.902 &| z + [ a4 + 0.1495 /3 4 ,o - 167.43] b BZ + 



(4.64) 



For Nf = 4 we obtain the following numerical values at the first three truncation orders: 
a ffs °) /A) = 0.312, 0.279 and 0.488, corresponding to a^ l s (/i 2 -»• 0) = 0.649, 

0.582 and 1.017, respectively. The higher-order corrections are certainly significant, in 
particular the NNLO ones. The results suggest a fixed point value that is marginally 
perturbative. This is consistent with the conclusion we reached in Sec. 14.31 by considering 
directly the fixed-point solution of the effective-charge beta function. 



5 A dispersive representation of the moment— space 
Sudakov exponent 

In Sec. 14. II we have written a dispersive representation for the physical Sudakov anomalous 
dimensions. Its ingredients have been studied in detail, the characteristic functions in 
Sec. [3]and the corresponding effective charges in Sec. HI We are now ready to use these tools 
to construct a dispersive representation of the Sudakov exponent based on the definitions 
in Sec. [21 We proceed as follows: in Sec. 15.11 we present a short derivation of the main 
result, a general dispersive representation of the Suakov exponent. Then, in Sec. 15.21 we 
present explicit results for the moment-space characteristic functions entering the Sudakov 
exponent in various examples and analyze their limits. Finally, in Sec. 15 .31 we return to the 
general discussion by considering in detail the role of virtual corrections in the large-TV 
limit. This discussion elucidates the interpretation of the dispersive formula. 



5.1 The Sudakov exponent: a general dispersive representation 

The crucial step here is going from momentum space, where the real-emission contribution 
is finite and can be considered separately, to moment space, where an infrared singularity 
is generated, which requires cancellation with virtual corrections. This cancellation is 
incorporated in the Sudakov exponent of Eq. (12. 7p . or equivalently in f)2.12p or (12.151) . 
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It is convenient to consider first the Sudakov evolution equation obtained through the 
logarithmic derivative of (12.151) : 
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(5.1) 



where in the second line we have inserted the derivative of the momentum-space kernel 
(13. 2p that admits infrared factorization. Owing to the cancellation between the real and 
virtual terms for r — > 0, the 'jet' and 'soft' contributions to (15.11) are separately infrared 
finite. This explains the third line. The purely real Laplace integral, 



f°° dr 

I «' —Nr rrt 2 

/ — e J[rm 

Jo r 



however, is infrared divergent. Here, the dispersive integral becomes handy. Substituting 
into (15.11) the dispersive representation of the Sudakov anomalous dimensions according 
to (14.11) and (14. 2p . respectively, we obtain: 



din Sud (m 2 , iV) C R 



dlnm 2 



Po 



00 dr 
r 



00 



dr 
r l 



e~ Nr - e~ r 



-Nr 

e — e 



de 



-pj{em 



de 



Ps(em 2 ) f s (e/r 2 )-f s (0) 



(5.2) 



where, as usual, e = fi 2 /m 2 . Next, changing the order of integration and defining the 
moment-space characteristic functions by the Laplace transform of the momentum-space 
ones: 



Qj (Ne) 



dr 



e-"'*>(e/r) 



Jo r 

we obtain an elegant dispersive representation of the Sudakov evolution equation: 



(5.3) 



dlnSud(m 2 ,iV) C R 



dlnm 2 



de 



/ - Pj (em 2 )[gj(eN)-gj(e) + lnN 



00 de 
e 



Ps(em 2 ) (g s (eN 2 ) - Q s (e) + In iV 



(5.4) 
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where the IniV terms originates in the terms proportional to J~j(0) and Fs(0) in (15.21) . 
where we have evaluated the trivial integral 



/ y[e-""-e-']=-l„JV 
Jo r 



and used Eq. (I3.27P setting ^j-(O) = ^s(O) = 1. We have explicitly written the 'jet' and 
'soft' contributions to the r.h.s of H5.4[) as two separate integrals, in order to emphasize 
the fact that they are separately finite: the combinations 

Qj (eN) - Qj (e) + In N and Q s (eN 2 ) - Q s (e) + In N (5.5) 

are each falling as a power of e; indeed the expansion of the characteristic functions at 
small arguments starts as 

Qj (eN) = - ln(eiV) + J x - lE + 0(eN) 

Qs (eN 2 ) = -~ ln(dV 2 ) + l - Sl - lE + ((eN 2 ) 1 / 2 ) , ^ 

where Si and Ji are the log-moments defined in (14.121) . These expansions are most easily 
derived using Eq. (16.191) below. Explicit examples will be given in Sec. 15.21 

Because 'jet' and 'soft' contributions to the r.h.s of (15.41) are separately finite, evolution 
equations similar to (15 .4p can be written separately for the 'jet' and 'soft' Sudakov factors 
(see e.g. Eqs. (8) and (22) in Ref. [54]) whose product is Sud(m 2 ,iV). The definition 
of such factors, however, requires to introduce a factorization scale. This will not be 
necessary in the derivation we present below, where we shall combine the two dispersive 
integrals in (15.41) into one. 

Note that integrating Eq. (15. 4p by parts we obtain an alternative expression in terms 
of the Sudakov effective charges: 



dlnrn 2 ft | J 6 V 

(5.7) 

- ar k (em 2 ) (q s (eN 2 ) - Q s (e) 
Jo e v 

The final step in deriving a dispersive representation of the exponent is to integrate 
the evolution equation (15.41) . Quite conveniently, the entire m 2 dependence of the r.h.s in 
(15 .4p is through pj(em 2 ) and ps(em 2 ). Using (14. 3p we can therefore readily integrate (15 .4p 
getting: 



Sud(m 2 , N) = exp -| — - / — 

Po Jo e 



' ^ -j a^(em 2 ) [Qj (eN) - Qj (e) + In n] 



aT k (em 2 ) (q s (eN 2 ) - Q s (e) + IniV 



(5i 
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which is the main result of this section. Explicit results for the moment-space character- 
istic functions Gj,s i n a few examples will be computed and analyzed in the next section; 
the final expressions are compiled in Table HI 

The exponent ( 15. 8p is finite owing to the exact cancellation of logarithmic singularities 
to any order in perturbation theory; there are two levels of cancellation: first between real 
and virtual terms within the separate 'jet' and 'soft' parts, and second between these terms. 
The latter cancellation mixes the infrared and ultraviolet in an interesting way. To see this 
note that for e — * the convergence of the separate 'jet' and 'soft' parts in (15. 8p follows 
from the convergence of ( 15.41) (the power-like falloff of the combinations in ( 15. 5ft above). 
On the other hand, in contrast with (15.41) . for e — > oo the integrals corresponding to the 
'jet' and the 'soft' parts in (I5.8P are not separately finite. For the Gj,s terms convergence 
for e — > oo is guaranteed by the Laplace weight in (15. 3p . The ultraviolet divergence is 
therefore entirely due to the IniV term, identifying its origin as the cusp singularity. In 
conclusion, the In iV terms are required for the separate 'jet' and 'soft' integrals to converge 
for e — > 0, but render each of them divergent for e — > oo. Nevertheless, the exponent as a 
whole is finite owing to the fact that a<f nk (em 2 ) and a^ ink (em 2 ) are equal at leading and 
at next-to-leading order (see Sec. (14.11) ) and therefore their difference behaves as al(fi 2 ) 
at large fi 2 . 



process 


moment-space characteristic function, Gj{£) 


jet function 
(e.g. DIS) 


(l 4 - 5 Ei(U) + (f - |) e- £ 


process 


moment-space characteristic function, Gs{v 2 ) 


B decay; HQ 
Fragmentation 


- + V - cos(z/) - sin(z/) ^Si(z/) - ^ 


)- 


— sin(z/) + cos(z/) Ci(z/) 


Drell-Yan ; 
99 -> Higgs 




e + e~ — > jets 
C parameter 
(c = C/6) 


— - lE -\n(u) - - 3 F 2 ^ ll2 


16 ) 


7TZ/ 3 / i 3 u 2 \ 

+ 192 3F2 ^,2 2 ,| 16 J 


e + e~ — > jets 

Thrust 

(t = l-T) 


Ei(l,z/) 



Table 4: Summary of results for the moment-space characteristic functions Gj{£, — Ne) 
and Gs(v 2 = N 2 e) of some inclusive distributions. 



Note that a somewhat different picture emerges in the approximation where Oj ink (em 2 ) = 
aj ink (em 2 ) = a"g k (em 2 ), which is valid to NLL accuracy (see (I4.39P ) or in the large-/?o 
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limit. In this case, the IniV terms cancel, and we obtain 



C f°° df 

Snd(m 2 ,iV)| NLL = exp<j-^ / - <C(em 2 ) 



(Qj (eN) - Qj (e)) - (q s (eiV 2 ) - & (e] 



(5.9) 



Here each of the terms is separately finite for e — > oo, but the e — > singularities cancel 
in two levels: first between real and virtual, and second between the 'jet' and 'soft' terms. 



5.2 Moment— space characteristic functions and their expansions 

Let us examine now some explicit examples. Having determined the functional form of the 
momentum-space characteristic functions (Table [2]) we can readily compute the Laplace 
integrals ( 15. 3j) to determine the moment-space ones. The final results are collected in 
Table H 

In what follows we give further details on the calculation of the characteristic functions 
and address the convergence of the dispersive integral in the exponent (15.81) considering 
the asymptotic behavior of these functions at small and large e. In the Drell-Yan case 
we extend the discussion and compare the dispersive formalism developed here with the 
joint resummation formalism or Ref. [17]. For Drell-Yan we also derive an alternative 
representation of the exponent in terms of a Euclidean integral [68]. Such a representation 
does not exist in the other examples. 

The jet function 

Let us begin by constructing an explicit expression for the moment-space characteristic 
function of the 'jet', Qj(N/j, 2 /m 2 ), based on the momentum-space one given in (13.441) . 
According to ( 15. 3ft we have: 



SAO 



dr 



dr 



dw 



w 



-Nr -p I ^ 



J 



m 2 r 



-Nr 



6(r > /i 2 /m 2 

11 II 

~ 2w ~ 2^ 



1 - 



1 /i 2 



2 m 2 r 



1 + |4) E i( U)+(f-f ) e- 



(5.10) 



where £ = Nfi 2 /m 2 and where the Exponential-integral function is defined as usual by 

"°° da 



-of 



(5-11) 
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Let us now examine the convergence of the e integral in (15.81) . First, for e — > oo there 
is convergence for each of the separate Gj{Q terms owing to the behavior of this function 
for large £: 



GAS) 



3 1 



2£ 



1 



(5.12) 



In the e — > limit there is convergence in ( 15. 8ft only owing to cancellations between 
the separate terms. We obtain 



Gj(eN) =(- lE - ln(iV) - ln(e) - 



lE + \n(N) + \n(e)jjNe + 0(e 2 ), 

(5.13) 

making the combination Q j (eN) — Qj (e) + In N power suppressed as required for the 
integral to converge. 

The Quark Distribution function (B decay) 

Let us consider now the 'soft' function associated with the Quark Distribution in an 
on-shell heavy quark. Using in (15.31) the explicit expression for the momentum-space 
characteristic function in (I3.44p we obtain : 



r 

drj 



-UT) 



o V 



+ 



1 + T] 2 (1 + T) 



2\2 



+ - e -^Ei(l, -iv) + e lv Ei(l, iu) \ + —[ e" w Ei(l, -iv) - e i!/ Ei(l, iv 



1 

2 + 



V 



— cos(v) — sin(^) 



— sin(v) + cos(v) 



Ci(v) 



(5.14) 



where we have defined v = N/j,/m = N^Je and changed variables to rf = r 2 /e. The 
Exponential-integral function is defined in (I5.1ip above, and 



Si(f) = J j sin(t) , 

Ci(i/) = 7s + In(i/) + / 

Jo 



(5.15) 



(cos(t) - 1) 



Let us examine the convergence of the e integral in (15.81) . For e — > oo we find 



(5.16) 



which guarantees convergence separately for each of the Gs QD terms. As expected, in the 
e — > limit there is convergence in (15.81) only owing to cancellations between the terms. 
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We obtain 

Gs QI) (Nh) = (-7* - ln(iV) + \-\ ln(e)) + ^Ne^ - \n\ + ^Nh^ + 0(e 2 ) 

(5.17) 

making the combination <?s QD (^ 2e ) — &s QD ( e ) + hi iV power suppressed. 
The Drell-Yan soft function 

Next, let us consider the 'soft' function associated with Drell-Yan production. Using in 
(15. 3p the expression for the momentum-space characteristic function (13.491) . we obtain: 

^ = ^e-^(^) = \[^=^ = W), (5-18) 

where we have defined ^ = Nfi/m = A^^/e and changed the integration variable from r to 
( = 4e/r 2 . Here Kq is the modified Bessel function of the second kind. 

To examine the convergence of the dispersive integral in the exponent ( 15.81) for e — > oo, 
let us expand (15.181) at large v = Ny/e. The result is: 

fc D >^ie-^( 1 -iI + (-L)). (5,9) 

Obviously, the integrals over the Qs DY terms converge for large e. 
Next, consider the e — ^ limit. Here the expansion yields 

g SDY (N 2 e) = ( - In (NVe) - lE ) + ( - In (Nyft) - 1e + l)N 2 e + O ((iV 2 e) 2 ) (5.20) 

making the combination Q Sl)Y (N 2 e) — Ss OY ( e ) + hiiV power suppressed at small e, as 
required for convergence. Note that only even powers of v appear in this expansion. 

Let us return now to the physical Drell-Yan kernel. The Laplace equivalent of (13.141) . 
which differs from it only by 0(1/N) terms for iV — > oo, is 

d\nG qq (N,Q 2 ^ 2 F ) ^lnSud DY (Q 2 ,AQ d\nH DY (Q 2 , 

— ^ — " — — + — — + ° {1/N) ' (5 ' 21) 
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with 

d In Sud DY (Q 2 , N) 



ding 2 
2C F 

~K 

2C F 

~K 

2C F 
~K 
2C F 



2C, 



-Nr 



9(r < 1)) S DY (Q 



2„2\ 



dr 



(e~ Nr - 6(r < 1)) 



dp 2 

dp 2 
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dp 2 



(P 2 



Mink 

,2 a SvY 



p 



Mink 

,2 a SvY 



(p 2 



(p 2 



d 



dlnQ 2 
d 

'dlnQ 2 



o 



dp 1 - 

p2 



Ps DY (p 2 



p 2 

Q 2 



p- 



Q 2 r 2 
+ 7s + ln(iV) 



•(0) 



N 2 p 2 



1 



2C f 



2 °5dy^ 



2 In 



Q 2 

at (2v/ivv/g 2 

l + y/l- Ap 2 /Q 2 



>Sdy 



Q 2 



(5.22) 



-% 2 /Q 2 <l/4)x 



ln(p 2 /Q 2 



AT 2yW/Q 2 - ln(VgV^) 



cilnQ 2 

where in the second line we used the dispersive representation of S DY from (14.21) . in the 
third (see derivation in Sec. 15.31) we performed the integration over r and substituted the 
Laplace transform of JF DY according to (I5.3P (and used J r DY (0) = 1), in the fourth we 
performed integration-by-parts and in the fifth we inserted the explicit expressions from 
(I5.18P and (13.481) . Finally, in the last line we have changed the upper limit of integration 
over the Kq term, which modifies the result by an exponentially small correction at large 
N, and neglected the first term in the square brackets, 2 In ((1 + y/l - 4p 2 /Q 2 )/2) , which 

generates 0(A 2 /Q 2 ) power-suppressed corrections. 

Note that the logarithmic derivative with respect to Q 2 cannot be pulled out of the 
p 2 integral in (15.221) . which would then diverge. This divergence is the usual collinear 
divergence associated with the evolution of the quark distribution - - in the Sudakov 
exponent of infrared and collinear safe distributions, Eq. ( I5.37p . it cancels with the jet 
subprocess. In order to obtain Sud DY (Q 2 , N, fi 2 F ) one would obviously need to introduce a 
factorization scale. 

It is interesting to compare our final result for the large-iV Drell-Yan kernel (I5.22p 
with the equivalent expression derived in the joint resummation formalism [17]. To this 
end consider the result for the Eikonal cross section, Eq. (48) in [17], where dimensional 
regularization (D = 4 — 2e) is used to regularize the collinear divergence in the exponent: 



a q g(N,bQ,e) = exp 




uj q q(k 2 ,k 2 + k 2 ± ,fj 2 ,e)5(k 2 + k\-p' 



e"^ K (ViV 2 p 2 /Q 2 ) - In ^QV^ 



(5.23) 
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where b is the Fourier conjugate to the transverse momentum of the system, and uj qq - 
is the "web" [71] corresponding to radiation from the annihilating lightlike quark and 
antiquark, defined at fixed transverse momentum. One immediately recognizes that upon 
considering the small b limit, corresponding to the "total" cross section (or to a case 
where the resummation of transverse momentum logarithms is unimportant compared to 
threshold logarithms, bm <C N), Eq. (15.231) reproduces Eq. (15.221) up tJ3 terms that are 
finite at iV — ► oo, provided one makes the following identification^: 



%P 2 



lim <>V^) 



(5.24) 



where 

C F a^ Y V,/i 2 ,5) 



0o 



P 



jdk 2 J ( t^u qq {k\k 2 + kl^\e)5(k 2 + kl- V 2 ). (5.25) 



This gives a direct diagrammatic interpretation to the Minkowskian effective charge in 
terms of webs. 

Eq. (15.221) is the standard dispersive representation of the kernel that involves the 
Minkowskian effective charge. In the Drell-Yan case, the exists an alternative Euclidean 
representation [68]. To derive it one defines a Euclidean characteristic function by the 
following dispersion relation: 



QSr 



dy 



(5.26) 



Substituting this expression for Qs DY i n the s °ft P ar ^ °f ( 15.41) (after integration by parts) 
one gets: 



d\nSud DY (Q 2 ,N) 2C f 



dlnQ 2 



0o 
2C F 



d/j, 



Mink 

2 a SuY 



dy 



A 4 



1 + 



'<Sdy 



(yN 2 ) 



Ss£ (y) 



00 



(5.27) 



where in the second line we used the relation between the Minkowskian and the Euclidean 
effective charges, Eq. (14.41) . 

According to (15.181) . the explicit result for G§^(y), defined in f!5.26j) . is [68,94] 



(5.28) 



where Jo is the Bessel function of the first kind. 

' Finite terms are generated by a different upper limit on the momentum integration, Q 2 in (|5 . 23[) 
versus Q 2 /i in (|5.22[) . which amounts to a different split of the virtual corrections between the Sudakov 
factor and the hard function that multiplies it. 

■fin the e — > limit, the result is independent on the dimensional regularization scale /j,. For example, in 
case of a single gluon emission (Eq. (45) in [17]), the r.h.s in (|5.25p is lim e _> Cp(a s /w)(4:ire~' yE [i 2 /p 2 ) e /p 2 . 
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Finally, note that when using the Euclidean representation (I5.27p . power-like infrared 
sensitivity (i.e. infrared renormalons) is directly reflected in the expansion of the Euclidean 
characteristic function near zero. This should be contrasted with the Minkowskian repre- 
sentation where only non-analytic terms in the expansion are associated with renormalons 
(see Sec. [7] below). In the case of Drell-Yan production, only even powers of NA/m ap- 
pear [52], as can be seen [68] from the fact that Jq is an even function of its argument. 
Equivalently in the Minkowskian representation the non-analytic terms in the expansion 
of K appear with even powers (see Eq. (15.201) ) — this has been pointed out in [17,95]. 



The soft functions in event-shape distributions 

Let us compute the moment-space characteristic functions of the thrust and the C pa- 
rameter. 

The case of the thrust is simple: J-'stiz) = 0(z < 1). Using (15.31) we obtain 



Q St (y 2 



1 ^e-"/^ 



Ei(l, v) 



(5.29) 



The case of the C parameter is slightly more complicated. Let us write the momentum- 
space characteristic function as 



[z < A) + 6(z > 4) 



1 



1 



and split the integral accordingly: 



dz 



dz 



dz 
I — e 



Ei(l,i//2)+ f —e 



vx/2 



Z 
1 



^Ifz 



— -1E- ln(i/) - — 3^2 



1,1 

3 3 o 
2' 2' 



X' 

16 



192 6 



1 2 

L i 2 

2 ? £ 



V 

16 



(5.30) 

Let us examine the convergence of the e integral in (I5.8p . For e — > oo we find exponential 
suppression for the thrust and power suppression for the C parameter: 



Qstiy 



2 " ~ e"" 



QsJy*)~-- + -I 



7/- 



(5.31) 



In either case each of the Q$ terms in (15. 8p converge. 
For e-^0we find 



&> 2 



QsXJ 



My 
ln(z/) - j E 



1 2 1 3 

IE + v V H V 

1E 4 18 



7T 



1 



-1/ 



7T 

192 



+■ 0(" 4 ) ; 

v 3 + 0(u A ) 



(5.32) 
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making the combinations Gs(N 2 e) — £/s(e) + In N power suppressed at small e, as required 
for convergence. 



5.3 The large— N limit, strict factorization and finite terms 

As discussed in Sec. [2} the Sudakov exponent in a given process is unique as far as the terms 
that diverge at large-iV are concerned, but it is subject to convention in what concerns 
finite terms for N — > oo as well as any 0(1/ N) corrections. This has been illustrated 
by introducing three definitions for the Sudakoiv factor in Eqs. (12.71) . (12.121) and (12.151) . 
Our purpose here is to identify, in the framework of the dispersive approach, the unique 
large- N limit of the Sudakov exponent, which in contrast with Eq. (15. 9ft above, involves 
two scales only: m 2 /N and m 2 /N 2 (such a strict scale separation is a must in an effective 
field theory framework). We will also show that a convention-dependent part that violates 
this strict scaling is in fact necessary to render the Sudakov exponent finite. 



The large— N limit: disentangling N— dependent and constant terms 



Finite terms for N — > oo are generated by purely virtual corrections, which are propor- 
tional to 5(r) in momentum space. Purely virtual corrections are partially accounted for 
within the exponent and partially in the hard function H(m 2 ) in (12.61) or (12.141) . This split 
depend on the convention adopted, except for the requirement that all infinities are as- 
signed to the exponent where they cancel against the integral over the real-emission part. 
Let us make this dependence parametric by generalizing Eq. (12.121) though a modification 
of the cutoff value: 



Sud(m 2 , N, r ) = exp ^ Cr j — 



-Nr 



0(r < To) R(m 2 , r) 



(5.33) 



with 



T N {m 2 ) = H(m 2 , r ) x Sud(m 2 , N, r ) + AT N (m 2 ). (5.34) 

For r = 1 the constant term 0(N°) in the exponent matches the original Mellin-transform 
definition (12771) . 

It is straightforward to repeat the derivation of the dispersive representation with this 
definition. Proceeding as in (15.11) and (I5.2p we obtain: 



d\nSud(m 2 ,N, r ) 



d\nm 2 



Cr 
Po Jo 
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r°° dr 


Q -Nr_ 


'o r 




de 


r°° dr 


- Q -Nr _ 
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Jo r 





00 de 
e 



1 dR(m 2 ,r) 

6(r<r ) —\ ^ 

J a in m z 

— 6(r < r ) x 
?J (e/r) - ^(0)) - Ps (em 2 ) (f 8 (e/r 2 ) - T s (0)) 
Pj {em 2 )Hj{N, e, r ) - Ps {em 2 )H s {N, e, r Q ) 



(5.35) 
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with 



Hj(N,e,r ) 



H s (N,e,r ) 



dr r 



o 



r 



e- Nr -9(r<r ) [T a (e/r) - *>(0) 



^(JVe) - ^(e/r ) + ln(r ) + Te + ln(iV) 

"°° dr 



-jV> 



(r < r ) 



.F 5 (e/r 2 ) - T S (0) 



(5.36) 



£s(iV 2 e) - -^(e/r 2 ) + ln(r ) + 7 £ + ln(iV) 



where in the last step in (I5.35P we performed the integration over r substituting Qj(Ne) 
and Qj(N 2 e) according to (15.31) and used (13.261) to write 



" -^(e/r) = Tj{e/r Q ) ; 



o 



y^(e/r 2 ) = ^(e/r 2 ). 



Note that 1~Lj,s in Eq. (I5.36I) are infrared finite owing to the subtraction of 8(r < r ) in 
the square brackets; they are ultraviolet finite because of the differences Tj (e/r) — Tj (0) 
and Ts ( e / r2 ) — Fs (0)? respectively. Finally, performing the lnm 2 integration in (I5.35P 
we obtain the result for the Sudakov factor: 



Sud(m 2 , N, r ) = exp 



Pi 



de 



o Jo 



a^(em 2 )Hj(N, e, r ) - a^(em 2 )H s (N, e, r ) 



(5.37) 



Explicit results for H,j t s(N,e,r ) can be readily obtained by substituting Gj,s and Tj^s 
of Tables 0] and [H respectively into Eq. (I5.36P . 

Let us now check the convergence of the integral over e in Eq. (I5.37I) . Consider first the 
e — > limit. Using the small-e expansion of Gj,s in (I5.6P and that of J-j t s in (I4.17P we 
observe that the combinations in Eq. (I5.36P are power-suppressed at e — ► 0. Next consider 
the e — > oo limit. Gj,s an d Tj % s are ( & t least) power suppressed at large e. The remaining 
terms ln(ro)+7_e+ln(A^) are common to the 'jet' and 'soft' part, and therefore the integrand 
is proportional to the difference of the two effective charges, Oj ink (em 2 ) — a^ ink (em 2 ), which 
behaves as a^(/i 2 ) at large /x 2 . Thus, the integral (I5.37P is well-defined. 

Let us now isolate the iV-dependent terms in the exponent. This can be elegantly done 
by considering the limit r — > oo. Expanding TCj,s(N, e, r ) at large r we obtain: 



HAN,e,r ^l)=AG ( ; +v \Ne) + O{e/r ) 
H s (N,e,r » 1) = (iV 2 e) + O(e/r 2 ) 



(5.38) 
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where we used the expansion of J-j : s at small arguments, Eq. (I4.17p . and defined 



Ag ( I +v) (Ne) 



AQ { ; +v \N 2 e) 



r°° dr 


- Q -Nr _ { 


~ Jo r 




= Gj(Ne) 


+ ln(iVe) 


r°° dr 


' Q -Nr _ l 


~ Jo r 





F s (e/r 2 )- ^s(O) 



(5.39) 



g s (N 2 e) + ^ln{N 2 e)-^S 1 + lE 



where the superscript (r + v) indicates that in contrast with Qj t s of (15.31) . which consists 
solely of real emission contributions, through the subtraction of 1 in the square brackets 
(I5.39P includes also virtual corrections. This guarantees convergence for r — > 0. Note that 
AQj +v \Ne) and AQ^ +v \N 2 e) converge for r —> oo owing to the subtraction of ^Fj(0) 
and ^s(O), respectively. Thus, similarly to TCj t s(N,e,r ) these objects are well-defined. 
Based on these definitions and (15.31) it is straightforward to recover the expressions in 
terms of Qj(Ne) and Gs(N 2 e), respectively. 



In contrast to Hj t s(N,e,ro) the functions AQj +v \Ne) and A^s^(iV 2 e) depend on a 
single argument. Using (15.381) we therefore find that the iV-dependence of the jet function 
in (15.371) involves only the jet mass scale m 2 /N while that of the soft function only the 
soft scale m 2 /N 2 . This is the infrared factorization property we have already encountered 
in Sec. [3] when working in momentum space. We also observe, however, that despite the 
fact that the integral in (I5.37P is well defined for any finite r , the r — > oo limit does not 
exist: by taking ro — > oo inside the e integral a divergence is generated for e — > oo. Thus, 
the O(e/r ) and 0(e/r 2 ) contributions that distinguish Hj,s(N,e,r ) from their simple 
ro — > oo limits (thus violating strict factorization) are in fact necessary for convergence. 

Having separated the iV-dependent terms from the cutoff-dependent constant terms, 
Eq. (I5.37P can be elegantly written as follows: 



C R r°°de 



Sud(m , N, r ) = exp <( — / — 



1 



- aT\em 2 ) ( AG^(N 2 e) - ^ s (e/r 2 ) 



where for the r dependent terms we used the following definitions 

r e / r o dv r • • i 

AFj{e/r ) = - - Tj{y)-Tj (0) = Tj (e/r ) + In (e/r ) - h 

Jo y L J 



t/4 



dy 

y 



Ts (y) - T s (0) = T s (e/r 2 ) + In (e/r 2 ) - S x . 



(5.40) 



(5.4i; 



Note that Ag ( j +v) (Ne) and Ag£ +v) (N 2 e) in (15^91 are essentially Gj(Ne) and Gs(N 2 e), 
respectively, minus the leading terms in their expansions, Eq. (15.61) . Similarly AjFj(e/r ) 
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and AJF^c/Yq) in (15.411) are Tj{ejr§) and JF^c/tq) minus the leading terms in their 
expansions, Eq. (I4.17p . Consequently, the convergence of (15.401) for e — > is guaranteed 
for each of the separate terms, which all fall as e 1 in this limit. In contrast, the convergence 
for e — > oo involves cancellation between all the terms. 



Re-derivation of the exponent using strictly factorized components 

We saw that upon removing the regulator r one recovers strictly factorized functions that 
fully capture the iV-dependent terms of the Sudakov exponent at large N. We can now 
repeat the derivation of the dispersive representation of the exponent with no constant 
terms, Eq. (I5.8p . using strictly factorized components. To this end let us define the formal 
'jet' and 'soft' Sudakov exponents: 

r°° dr 

Ej(m 2 /N) = C R - [e~ Nr - l] J(rm 2 ) (5.42) 
Jo r 

and 

E s (m 2 /N 2 ) = C R T - [e~ Nr - l] S{r 2 m 2 ) , (5.43) 
Jo r 

which include both real and virtual corrections. Owing to the virtual term (—1) these 
expressions are infrared finite, but since the integration extends to infinity they are ultra- 
violet divergent. In Appendix Owe show that these formal definitions correspond to taking 
the large iV-limit of the finite-iV physical evolution kernel, where the limit is taken with 
fixed Nfj, 2 /m 2 for the 'jet' and fixed N 2 p 2 /m 2 for the soft function. The characteristic 
functions of Eq. (15.391) are then identified with the corresponding large-iV limits of the 
full finite- N characteristic functions computed in Ref. [49]. 

The Sudakov evolution equation, Eq. (15. ip . can be expressed in terms of finite combi- 
nations of these quantities, where the ultraviolet divergences cancel out. We have 



dlnSudK,A0 = _ _ [Es{m 2 /N 2 ) _ Es{m2)] (5 44) 

Substituting into ( 15.421) and (15.431) the dispersive representation of the Sudakov anoma- 
lous dimensions according to ( 14. ip and ( 14.21) . respectively, we obtain the moment-space 
representations: 

Ej(m 2 /N) = ^ * Pj (em 2 )Ag^ (eN) 

P ° J ° £ (5.45) 

E s (m 2 /N 2 ) = ^ / - Ps (em 2 )Ag^ (eN 2 ) , 
Po Jo e 

where, as usual, e = fi 2 /m 2 , and AQj +v ^ (eN) and AQ^ +V ^ (eN 2 ) have been defined in 
( I5.39p . Note that here the integrands are well-defined and integrable near e — > but, as 
expected based on the definitions (I5.42p . the integrals diverge for e — > oo. 
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Integrating (15.451) by parts we obtain an alternative form in terms of the Sudakov 
effective charges: 



Ej(m 2 /N) = ^ I * a^(em 2 )AQ { ; +v) (eN) 



c 



E s (m 2 /N 2 ) = ^ f° - ar\em*)Agt v) {eN 2 ) . 
Po Jo e 



(5.46) 



Next, substituting in Eq. ( 15.441) the dispersive representations ( 15.451) of the jet and 
soft Sudakov exponents, we get the dispersive representation of the Sudakov evolution 
equation: 



d In Sud (m 2 ,iV) C R 



dlnm 2 



de 



de 

Mr+v) 



p s {em 



A^ +v) (eN) - A^ +v) (e) 

ag { ; +v) (eN 2 ) - ag [ ; +v) ( e 



(5.47) 



Using the expressions for AOj^' in Eq. (15.391) this yields Eq. (15 .4p . Finally, integration 
over lnm 2 gives: 



Sud(m 2 ,iV) = exp{^ [°°- 



a M j*\em 2 ) I A^ +v) (eN) - Ag^ v> (e) 



a^ k (em 2 ) ( Ag£ +v) (eN 2 ) - A^ +v) (e) 



(5.48) 



which yields Eq. (15. 8ft . Thus, we have recovered the final result of Sec. (15. lft in terms 
of the characteristic functions Ag'j^ that include virtual corrections. This derivation 
elucidates the origin of the ln(iV) terms in Eq. (15.81) . Eq. ( 15.481) also exhibits explicitly the 
fact that the integral converges for e — >• 0: here each term is separately power suppressed 
in this limit. Convergence in the e — > oo limit relies on cancellation between the terms. 



6 Relation with the Borel representation of the Su- 
dakov exponent 

The Borel formulation of Sudakov resummation has been the subject of much theoretical 
work over the past few years [5,22,23,32,43,54]. It led to a successful phenomenology in a 
range of applications. The purpose of the present section is to study in depth the relation 
between this formulation and the dispersive one, which was presented in the previous 
sections. We will show that the two provide different ways of summing up the same set 
of radiative corrections. In the next section we shall address potential differences at the 
power level. 
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6.1 Borel representation of the anomalous dimensions 

Let us begin by introducing the scheme-invariant Borel representations of the Euclidean 
effective charges (14.51) : 

r°° ( a 2 \ u 

arV)= / duT{u) — B[a^]{u) 

J ° V/i 7 (6.1) 



r°° /a 2 \ u 

c#V) = J duT{u) B[ay] 



where A 2 is defined in MS. We emphasize that the Borel integrals do not necessarily 
converge for arbitrary /i 2 : all that is required is that they exist for sufficiently large fi 2 . As 
usual, while the integrand in (16.11) is scheme invariant [64], the separation between T(u) 
and the Borel function B [a^] (u) itself is scheme dependent. In practice (see e.g. Eq. (2.8) 
in [41]) it is convenient to use the 't Hooft scheme where T(u) absorbs the dependence on 
the two-loop coupling 5 = (3i/f3 2 , 

1 (1 + uo) 

Note that in the large-/?o limit T(u)\ lSiTgef3o = 1 and 

B [«T] («)L A = B[aT] (u)L ge , o = el" (6.3) 

According to (14. 4p the Borel representations of the timelike couplings corresponding 
to (16.11) are obtained by analytic continuation: 



/*oo / A 2 \ u 

Jo ™ 
a^V) = [duT{u) (£f ^B[ar] («)• 



(6.4) 



Substituting these expressions into (14.11) and (14.21) . respectively, we obtain the Borel rep- 
resentation of the physical Sudakov anomalous dimensions of ( 13. 31) [5,54]: 



1 r°° /A 2 \ u 

1 poo / \ 2 \ U 

where 

S<7 («) = S [a^ 1 ] («) x Mj(u) 



(6.5) 



Sin 7TM 
7TU 

sin 7tm 



(6.6) 



S 5 («) = - -B[af«](u) x B 5 ( M) 

7TM 
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with 



lo z Jo z U kl 



(6.7) 



where Fj,s{.z) are the characteristic functions defined in the large-/?o limit ( 13. 23ft and 
Jfc and 5^ are the corresponding log-moments defined in (14.121) . Thus, the expansion of 
Mj t s(u) also provides a convenient way to extract the log- moments. Obviously, the same 
information about higher orders and infrared sensitivity (power corrections) contained in 
the characteristic functions J-j } s is encapsulated by the corresponding Borel functions 
Mj s(u). Table [2] summarizes the explicit results for Mj t s(u) in various inclusive distribu- 
tions. 

Note that 

^ / \ i sin itu 5 . , 

BjAu)L rge(3o = ea" x B JiS ( u ). (6.8) 

Although Mj t s(u) have singularities at certain integer values of the Borel variable, these 
do not generate renormalons at large /3q in the anomalous dimensions J{k 2 ) and S(k 2 ) 
because of the sin(7ra) factor associated with the timelike nature of these quantities. 

We further expect that B [a^s] (u) and thus, in virtue of (16.61) . also Bj^s{u) are 
renormalon free at finite (3q. Absence of renormalons in anomalous dimensions is expected 
on general grounds^, but a formal proof in the non-Abelian case is still absent and would 
be very important. 

When using the Borel formulation, the deviation of Bj t s{u) from their large-/?o limits 
is computed order-by-order in perturbation theory. It amounts to a multiplicative modi- 
fication of the large-/?o result [5] by a factor Vj t s(u), see e.g. Eq. (19) and (26) in [54], 
which we can now readily identify as 

VjAu) = J^yJ =B[a2$W (6-9) 

^[ a J,s\ W I l ar g e ft, 

This gives a new perspective on the physical meaning of this function. 



6.2 The Sudakov exponent 

To demonstrate the correspondence between the dispersive representation of the Sudakov 
exponent and the Borel one, let us now re-derive the Borel formula for the exponent 
(cf. Eq. (27) in [54]) following the steps of Sec. E] above. We will then present a second 
derivation starting from the dispersive formula (15.81) . 

§Note however we deal here with physical anomalous dimensions, for which this expectation might not 
be realized. 
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Direct derivation 



Inserting the Borel representation of the Sudakov anomalous dimensions in ( 16. 51) into ( 15. II) 
we obtain: 

dlnSud(m 2 ,AO 



dlnm 2 



c R i * 



J(/-/,r')-C K / £ 



-TVr 

e — e 



Cr 

A) I Jo 





e _7v r _ e _ r 






r 




<ir 


e _jVr _ e -r 






r 





T(u) 



duT(u) 



A 



2 \ " 



A 2 \" 



r 2 m 2 



Bs(u] 



(6.10) 

where the two terms corresponding to the 'jet' and 'soft' contributions to the r.h.s are 
separately finite. Here one can change the order of integration using u as an infrared 
regulator. Using 

"°° dr 



-Nr — u 



N u T(-u) 



we obtain an elegant Borel representation of the evolution equation: 



dlnSud(m 2 ,iV) C R 



dlnm 2 



c { r°° ( a 2 \ " 

^|jf duT(u)(— 2 j Bj(u)T(-u)(N u -l) 



(6.11) 



duT(u) 



A 5 



5 5 (?z)r(-2u)(Ar u - 1 



where, as before, the two Borel integrals are separately finite, for u — > 0. We note that the 
terms proportional to N u and A^ 2n in Eq. (16 . 1 1 j) . respectively, are the Borel representations 
of the formal 'jet' and 'soft' Sudakov exponents (I5.42j) and (15.43p . where u serves as a 
regulator. 

It is straightforward to see the relation with Eq. ()5.4|l N: 



- Pj (em 2 ) ( Qj (dV) - Qj (e) + In N 



-a- nk (em 2 )(^(eiV)-^( e ; 



A 



2 \ u 



duT(u) — Bj(u)T(-u)(N u -l) 



(6.12) 



and 



j Ps(em 2 ) (p s (eN 2 ) - Qs (e) + In JV 



°° * ar(em 2 ) (g s {eN 2 ) - £ 5 (e) 



duT(u) 



A 



77T 



2 \ « 



£ 5 (u)r(-2 W )(Ar 2 "-l). 



(6.13) 



^Note that upon using (16. 4p . Eqs. (|6.12p and (|6.13p can be derived from (|6.16|) below. 
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Note that in (16. lip . (16.121) and (I6.13P one must assume some prescription for the Borel 
singularities at positive integer and half integer values of u. 

The final step in the derivation of the Sudakov factor is to integrate (I6.1ip over 
lnm 2 . As in ( 15.81) this requires to combine the two integrals over u in ( 16.111) . The re- 
sult (c/. Eq. (27) in [54]) is: 



Sud(m , N) = exp < — 



C R f°° du 
Po Jo u 



/ A 2\« 

T{u) — r 



Bj{u)T{-u){N u -l) 
B s {u)T{-2u){N 2u -l) 



(6.14) 



This is the Borel equivalent of ( 15. 81) . 
Derivation based on the dispersive formula 

An alternative route is to insert the Borel representation of the Minkowskian couplings 
(16 .4p directly into Eq. (15.81) . This yields: 

C R f°°de 



Sud(m , N) = exp 



A) Jo e 

A 2 \« 



duT(u) 



sin tiu 



7VU 



poo 




i duT{u) 


( A2 ) 


lo 


\em 2 J 



sin ttu 



B[ay l ] (u) [g J {eN)-g J {e)+hiN 
B [af d ] (u) (g s (eiV 2 ) - Q s (e) + In N 



(6.15) 

Upon combining the two Borel integrals and changing the order of integration the lniV 
terms cancel and the e integrals over the Gj,s terms take the form: 



00 de 

-e- u [Gj(eN)-Qj(e) 



o 



-(N u - l)T{-u)Mj(u) 



u 



r -e- u (g s {eN 2 ) - Q s (e) ) = --(N 2u - l)T(-2u) B s (u) 
J t V / It 



(6.16) 



where we have computed the integrals using the relations (15. 3p and (16. 7p . which together 
yield a direct relation between the moment-space characteristic functions Qj t s ( e ) an d the 
Borel functions Mj t s(u) (all functions are defined in the large-/?o limit): 



Jo y 



'o y 



— y Fj{y) 
o y 



°° de 1 

-e~ u gj(e) = —T(-u)®j(u) (6.17) 

o e u 
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and 



-u 



z 



- t- u Q s (e) = --Y{-2u) M s (u) (6.18) 
e M 



It is straightforward to check this correspondence in the explicit examples considered 
above. Using (16.161) in (16.151) we immediately recover the known form of the exponent in 

(EH- 

Finally, it is useful to write the inverse relations to (I6.17P and (16.181) : 



Qs(e) - 



27T2 



2m 



30 dv 

— e u T(-u)®j(u) 

oo U 

— e u T(-2u) M s (u). 
u 



(6.19) 



Using these relations with the u — > expansion of Mj^(u) in (16. 7p we can easily derive 
the O(e ) term in the expansion of Gj,s{ e ) i n (I5.6P by applying the Cauchy theorem. 



7 Power corrections in the Sudakov exponent 

As discussed in the introduction infrared renormalons are present in the Sudakov exponent. 
This is evident in Eq. (I6.14p . As always, infrared renormalons reflect sensitivity to large- 
distance dynamics and indicate the presence of non-perturbative power corrections. From 
(I6.14p one can immediately deduce that these power corrections (1) exponentiate together 
with the logarithms; and (2) are enhanced at large N: the 'jet' and 'soft' contributions scale 
as integer powers of NA 2 /m 2 and NA/m, respectively. Such parametrically-enhanced 
power corrections appear exclusively through the Sudakov factor and have a significant 
impact on the distribution in the threshold region. 

The renormalon technique offers here a unique window into the non-perturbative side 
of the problem. Renormalon analysis allows us to find the specific pattern of power 
corrections [54]: in a given process only certain power corrections appear while others 
are absent. This pattern of power corrections is linked with the symmetry properties of 
the source of soft gluon radiation, and it is very much process dependent. The absence of 
specific power terms can be seen as an extension of the notion of infrared finiteness beyond 
the logarithmic level, i.e. "infrared safety at the power level". 

In the Borel formulation the presence of infrared renormalon ambiguities is transparent: 
the Borel integral in (I6.14p is obstructed by simple poles at integer and half-integer values 
of u, except where the Borel transforms of the Sudakov anomalous dimensions Bs{u) and 
Bj(u) of (16.61) vanish. It is natural then [22,54] to take the Principal Value of the integral 
in (I6.14p as a definition of the perturbative sum, and use the renormalon ambiguities as a 
basis for parametrization of power corrections. 
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In the dispersive approach, power-like sensitivity is somewhat less obvious. Consider 
for example the soft (S) Sudakov factor as can be obtained from (15.41) through integration- 
by-parts (an example is provided by Eq. (15.221) for the Drell-Yan case): 



dlnSud(m 2 ,iV) 



dlnm 2 



C 



de 



^ / -cCW) Gs{N 2 e)-g s (e) 



[m 2 ,N). 



(7.1) 

Only non- analytic terms in the small-e behavior of the Minkowskian characteristic function 
Qs (e) are related with renormalons [45,49,58], indicating corresponding power corrections. 
In a Euclidean formulation, 



dlnSud(m 2 ,iV) 



dlnm? 



Cr 

A) Jo 



dy 



ar\ym 2 ) Gr 1 {N 2 y) - Gr 1 (y) 



(7.2) 



such as the one of Eq. (I5.27P in the Drell-Yan case, the relation is simpler: all the terms 
in the small-y expansion of a^ ucl are related with renormalons. Although a Euclidean 
formulation (17. 2p does not exist in general — amongst the examples considered in this 
paper it exists only in the Drell-Yan case — it is instructive to consider it first since, 
in contrast with the Minkowskian one, it facilitates separating between the infrared and 
ultraviolet contributions in a straightforward manner by splitting the momentum integral: 



dlnSud(m 2 ,iV) 



din m 2 



< cl (^ 2 , N- tf) + Ei: cl {m\ N- /i 2 ) ee £ Eucl (m 2 , N) 



where 



E- cl (m 2 , N; /i 2 ) = I ^ ar\y ™ 2 ) {N 2 y) - Qf* (y) 



dy 

Po Jo y 



(7.3) 



(7.4) 



where A <C hi <C Q/N is a momentum cutoff and yi = \x\jm 2 . Using (17.21) and (17.31) the 
remainder is given by 



^ I ^ar\ym 2 ) Qf* (N 2 y) - fig- (y) 



dy 



(7.5) 



One can now expand the Euclidean characteristic function at small y under the integral 
(I7.4p . getting the power corrections. According to Eq. (I5.26P the Euclidean characteristic 
function (?J ucl (y) is the integral over the discontinuity of the Minkowskian one: 



dQT\y) 

dhxl/y 



— lm\g s (-y-iO) 

7T 



/->Eucl 



, °° del 

(y) = - I --lm\g s (-e-tO) 



(7.6) 



This is why only non-analytic terms in the expansion of the Minkowskian characteristic 
function Q$ ( e ) are relevant. 

The extension to the general case, where (17.21) does not exists, is based on ( 17.3ft where 
E^ cl (m 2 , N; /i 2 ) is still defined by (17.4ft but E^ cl (m 2 , N; n 2 ) does not have a representation 
of the form (f7.5p ; it must instead be defined by first reverting to a Borel representation 
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[45,58]. Since in this section we are focussing on power corrections, we do not write down 
the explicit form of the latter, but instead refer the reader to [45,58]. There is one point 
that needs to be stressed, however: Eq. ( 17.31) will not coincide with the Minkowskian 
dispersive integral (17.11) unless the Euclidean Sudakov coupling is causal. 

A causal coupling admits the dispersion relation (14.41) : it is analytic except for the 
cut on the timelike axis. In this case the entire first sheet of the complex momentum 
plane is mapped onto a compact domain in the complex coupling plane whose boundary is 
the image of the time-like axis (the simplest example where this is realized is a two-loop 
coupling with negative 0i/Po, see Fig. 2 in Ref. [74]). The three- and four-loop coefficients 
computed in Sec. (14. 31) in a few examples, have been found to have opposite signs to the 
one- and two-loop ones, opening the possibility that a non-trivial infrared fixed-point 
appears in these effective charges within perturbation theory. While the stability of this 
fixed point and the presence or absence of Landau singularities obviously depends on higher 
orders, a causal structure becomes possible. In this scenario the dispersion representation 
(I4.4p holds for the perturbative coupling itself, which does not need to be modified by any 
power terms, and then 

E Mink (m 2 ,N) = E Encl (m 2 ,N). 

E Mink (m 2 , N) still differs from the Principal Value of the Borel sum. But this difference is 
then entirely due to power corrections that are related to infrared renormalons [58,97-99], 
and genuinely reflect sensitivity to long-distance dynamics. The physical distribution is 
expected to differ from this perturbative result by power corrections whose pattern follows 
the renormalons. 

In contrast if a§ ucl (/c 2 ) has Landau singularities in the first sheet of the complex k 2 
plane, as indeed occurs in the large-/?o limit (see Appendix [Bj, then 

E Mink (m 2 ,N) E Eucl (m 2 ,N), 

and, although E Mink (m 2 , N) in (17. ip is still finite and yields a unique result, the dispersion 
integral (17.11) differs from the corresponding Borel sum by power corrections that are not 
related with infrared renormalons [45,49,58,61,62]. 

In the presence of Landau singularities in the Euclidean couplings, the Minkowskian 
couplings cij ll g(k 2 ) are still finite in the infrared limit. One may therefore be tempted 
to consider the result of (15.81) or even the approximation of (15.91) as a "non-perturbative 
model"|l]. One should be aware, however, that this result involves additional power cor- 
rections that are not associated with long-distance dynamics. As we recall in Appendix 
|B| using the dispersive integral amounts to "analytization" of the coupling that removes 
the Landau singularity. This is why the dispersive integral differs from the Borel sum by 
power terms that are not all related to renormalons. In the large-/?o limit this difference 
is given by the second term in the curly brackets in (IB. 101) . An explicit example is given 
in Appendix [D] 

In the following we derive the small-e expansion of the Minkowkian characteristic 
function, and then deduce from it the corresponding expansion of its discontinuity, namely 
of the Euclidean characteristic function. 

related model has been proposed in [39], where the fixed-logarithmic-accuracy formula for the 
exponent has been rewritten in terms of the analytic coupling (|B.8|) . 
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The relations (16.191) can readily be used to see that the singularities of T(—u)Mj(u) 
and r(— 2u) Ms(u) at integer or half integer values of u, u = k/2 (k is a positive integer), 
are associated with terms of order 0(e k ^ 2 ) in the expansion of the characteristic functions 
Gj,s(t) a t small e. Non-analytic terms scaling as e fc//2 in this expansion are related to 
infrared renormalons at u = k/2 and therefore indicate 0((A/m) k ) power corrections in 
the exponent. Owing to the iV-dependence of the exponent (15. 8p through Qj(Ne) and 
Q s (N 2 e), power corrections associated with the 'jet' and 'soft' subprocesses will scale as 
(Ne) k / 2 and (N 2 e) k ^ 2 , respectively. 

To see the precise relation between Borel singularities in (I6.14p and non-analytic terms 
in ( 15.81) . let us first recall the Borel singularity structure and then translate it into the 
terminology of the dispersive approach. The universal T(—u) and T(—2u) factors give 
rise to poles in (j6.14j) . The functions Bj t s(u) multiplying them may nevertheless vanish, 
cancelling the potential pole. According to (16. 6p this depends on the properties of Mj^ s (u) 
as well as the Borel transform of the corresponding Sudakov effective charge B \alj 5] (u). 
We assume, as in previous work [54], that the latter is innocuous: in the large-/3 limit 
it is given by (16. 3p . so it has no poles nor zeros at integer or half-integer locations, and 
we assume that singularities or zeros will not develop there beyond this limit. In contrast 
Mj^s(u) may have poles at integer values of u, as well as zeros at integer and half integer 
values of u. Note that the assumption above implies that the zeros of Bj S {u) present at 
large-/?o, which cancel some potential renormalon singularities in the exponent, are still 
present in the full theory. 

Thus, according to (16.191) non-analytic terms in the characteristic function emerge in 



the following cases: 

• In when Mj(u) has a pole, i.e. when T(—u)Mj(u) has a double pole. This 

occurs at u — 1 and u = 2 (see Table [2]). Indeed, upon expanding the explicit result 
for Qj(e) we get: 



Me) ^ (-J E - ln(e) - ~) + ( 2 - i 7 - i He)) e + f - | + \ie + ~ hi(e) 



2 



(7.7) 



The non-analytic terms e hi(e) and e 2 ln(e) in this expansion indicates, respectively, 
NA 2 /m 2 and iV 2 A 4 /m 4 power corrections in the exponent (15.81) : see e.g. [5,57,68,96]. 

• In {?.s(e), at all half integer values of u, u = k/2 where k is an odd number, except 
when B^ii) vanishes. In this case Gs(t) develops a square-root singularity {\fe) k '. 
Thus, Gs{N 2 e) in (15. 8p gives rise to 0((NA/m) k ) corrections. 

• In £As(e), at integer values of u (u = k/2 with even k) where B^w) has a pole 
(so r(— 2u)M s (u) has a double pole). In this case Qs( e ) develops a logarithmic 
singularity: e fc//2 ln(e), leading to 0((NA/m) k ) corrections in (15. 8p . 

The three cases above are precisely the ones where -B.7,5(^)| lar g e ^ in Eq- (I6.8P does not 
vanish, i.e. where the singularities in T(—u) and T(—2u) are not cancelled. 
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Altogether the soft characteristic function has the following generic structure at small e: 

g s (e) = ~ ln(e) + \s x - lE + E Ck W e t + E Ck ^ ln ( e ) + anal y tic terms • 

(7.8) 



k odd 



k even 



or 



I k k k 

Gs{e) = - + 22 - C k (y/e) + 22 o °k £ k/2 ln ( e ) + analytic terms . 



2 2 ^2 

k odd k even 

Thus, upon extracting the discontinuity according to (17.6j) one obtains: 

k sin(/c7r/2 
2 

' " feodd 

or 



7] ^ = --Im^ 5 (-,) 



7T 



fc even 



E?^ fe/2 



(7.9) 



(7.10) 



f_l)(fc-l)/2 

sra/HE^— — (Vv) k 



E ^/ /2 = E c ^ fc 



/2 



(7.11) 

fcodd fceven fc=l 

The coefficients C*. of the first few non-analytic terms in specific examples are summarized 
in Table [5j The table clearly demonstrates that the renormalon singularity structure is 
process dependent. 



process 




£rn(£) 


e 2 ln(0 


e 3 ln(0 


t'HO 


e ho 


e HO 


jet function 
(e.g. DIS) 


J 


-1/2 


1/4 
















process 
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v 2 ln(z/) 




z/ 4 ln(i/) 




z/ 6 ln(i/) 


B decay; HQ 
Fragmentation 




It /A 





tt/24 


1/24 


-tt/160 


-1/360 


Drell-Yan ; 
99 -» Higgs 







-1 





-1/4 





-1/36 


e + e~ — > jets 
C parameter 


Sc 


tt/4 





tt/192 





tt/20480 





e + e _ — > jets 
Thrust 


s t 


1 





1/18 





1/600 






Table 5: Summary of results for the first few coefficients of non-analytic terms in the small- 
e expansion of the moment-space characteristic functions Qj{S, = Ne) and Qsiy 1 = N 2 e) 
of some inclusive distributions (cf. Tables H] and [2]) . 



The small-momentum expansion of the Euclidean characteristic function (17.111) can be 
readily used inside the integral in Eq. (I7.4p to obtain: 

E^(m\N^]) ^ -2jL f> x (N k - 1) x r^y^ariym 2 ) ■ (7.12) 

po k=l Jo y 



"V - 
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Obviously, in the infrared region the coupling may differ from its perturbative part. This 
is encapsulated in the moments. In this way the moment = O ((////to) fc ) can be used 
to parametrize non-perturbative power corrections of order O ((NA/m) k ) with any k in 
the Sudakov exponent. Recall that the k = 1 moment (with a choice of [ii = 1 or 2 GeV) 
has been used in [49] and following work (see e.g. [26-29]) to parametrize power corrections 
to event-shape distributions. The framework presented here allows one to identify this 
coupling to higher orders in perturbation theory as well as to parametrize in a similar 
manner the corrections that scale as O ((NA/m) k> ) using the higher moments. 

Assuming that the couplings have a causal analyticity structure one can directly 
parametrize these functions in the infrared consistently with the dispersion relation (14.4R as 
well as their known ultraviolet evolution of Eqs. (14.361) and (14. 37ft ; a simple one- loop model 
of this kind has been considered in Ref. [100]. Upon taking the time-like discontinuity to 
obtain a^ in g one can readily evaluate the Minkowskian integral in the exponent, Eq. (15.81) . 
getting the all-order sum already including power corrections of the renormalon type. An 
even simpler possibility would be to parametrize directly the Minkowskian couplings a^ 1T g 
in the infrared, consistently with their ultraviolet evolution^]. In this way one can set up a 
simple yet fully consistent "power-correction phenomenology" without ever dealing with 
individual powers and without introducing any momentum cutoff. 

8 Conclusions 

We have presented a general formalism for Sudakov resummation based on dispersion 
integrals. The expression for the Sudakov factor in infrared and collinear safe distributions 
involving a jet function (J7) and a soft function (S) is summarized by Eqs. (15.480 and 
(15. 8p or (I5.40p . This formulation consists of two ingredients [68,69]: the first, Qj^s-, are 
characteristic functions that are computed analytically in the large-/9o limit and encode 
information on power corrections; the second, aj^, are Sudakov effective charges that are 
defined order-by-order in perturbation theory and encapsulate the non-Abelian nature 
of the interaction. In what concerns the soft function (S) both these ingredients are 
process-dependent owing to the fact that large-angle soft gluon radiation depends on the 
space-time geometry of the hard partons that radiate. In contrast, there is a unique jet 
function (j7) characterizing soft and collinear radiation from a jet with a fixed invariant 
mass. Explicit results for the characteristic functions in a variety of processes have been 
compiled in Table IU and the corresponding Sudakov effective charges have been computed 
and analyzed in Sec. HI 

The Sudakov effective charges are directly related to the physical anomalous dimen- 
sions. They can be computed to any order based on the conventional anomalous di- 
mensions defined by dimensional regularization. Furthermore, comparison with the joint 
resummation formalism [17] leads to a direct diagrammatic interpretation of the Sudakov 
effective charges in terms of 'webs', see (15.251) above. 

The dispersive approach presented here provides a realization of DGE [22,23,54,57]: it 

tin this case one should make sure that the resulting Euclidean coupling is causal at the perturbative 
level, and does not involve extra unwelcome power corrections, as in the APT example! 
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goes beyond the resummation of Sudakov logarithms per se by incorporating an internal 
all-orders resummation of running-coupling corrections. This resummation guarantees 
renormalizat ion-group invariance. Owing to the enhancement of subleading logarithms 
that are associated with the running of the coupling [22,57], this additional resummation 
leads to a significant improvement over the conventional approach to Sudakov resumma- 
tion where a renormalization-scheme-dependent truncation is performed, guided solely by 
the logarithmic accuracy criterion. Beyond the perturbative level, the present approach 
facilitates a systematic analysis of power corrections based on renormalon ambiguities 
which reveal themselves through the discontinuities of the characteristic function Gj,s- 

We have shown that there is a direct correspondence between the scheme-invariant 
Borel formulation and the dispersive one, and derived all-order relations between their in- 
gredients. The two formulations capture the same set of radiative corrections. Yet, at the 
power level they provide different regularizations of the sum. As far as renormalon singu- 
larities are concerned, this difference does not pose any difficulty: the two regularizations 
are in principal equivalent, as they differ by power corrections of the same parametric 
form that need be introduced to account for genuine non-perturbative effects. In con- 
trast, if the Sudakov effective charges have Landau singularities, as occurs for example in 
the large-/?o limit, the dispersion relation is violated. Then, the dispersion integral differs 
from the Borel sum by additional power corrections that are not related with renormalons. 
This complication does not arise if the Sudakov effective charges have a causal analyticity 
structure, making them consistent with the dispersion relation. 

The dispersive approach offers a convenient way to parametrize the non-perturbative 
power corrections exposed by the renormalons by means of integrals over the coupling in 
the infrared region [49]. We find that in the context of Sudakov resummation, the infrared- 
finite-coupling approach by Dokshitzer, Marchesini and Webber is of special interest: in 
contrast with the general situation here the coupling — the Sudakov effective charge - 
can be systematically identified to any order. Moreover, as discussed in Sec. I4.3[ there 
are indications that these effective charges may reach a finite limit in the infrared already 
within perturbation theory. Remarkably, despite the fact that the evolution of the Sudakov 
effective charges becomes process-dependent at three loops, the infrared limit itself turns 
out to be universal: it depends only on the cusp anomalous dimension. 

A particularly attractive example from the point of view of the dispersive approach is 
that of Drell-Yan or Higgs production. This has two aspects: first, there exists a Euclidean 
representation (I5.27P where the identification of large-distance effects is transparent, and 
second the relevant Sudakov effective charge admits an infrared fixed point already at the 
perturbative level as both the three-loop and four-loop coefficients of the effective-charge 
beta function are negative, see Eq. f)4.58p . 

A significant effort has been put in the past decade in examining the hypothesis of a 
universal [26-28,49] infrared-finite coupling in the context of event-shape distributions. 
Experimental data were primarily used to test the assumed universality of the first power 
moment of this "effective coupling". Indeed, overall this assumption is supported by 
data [29]. Despite the different resummation formalism utilized in these studies compared 
to the one presented here, the notion of the infrared-finite coupling, and the way it re- 
lates to power corrections, are the same. So far this "effective coupling" has only been 
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identified to NLO, and its universality has not been established theoretically. This has 
now changed: the present formalism uniquely identifies the effective charges relevant for 
Sudakov resummation to all orders in perturbation theory, and shows that their univer- 
sality does not extend beyond the NLO. It nevertheless shows that their infrared limit 
is universal. Approximate universality of power corrections, determined by the first few 
power moments of the coupling, is expected as a by-product. 

DGE, in its Borel formulation, has been successfully applied to phenomenology in a 
range of processes [54]. Most importantly, it extends the range of applicability of re- 
summed perturbation theory closer to threshold, far beyond what can be achieved by 
a conventional, fixed-logarithmic-accuracy approach. This should be attributed to two 
main factors: first the additional resummation performed, and second the possibility to 
identify correctly the pattern of power corrections, especially the first one or two renor- 
malon singularities. The main challenge has been to find an effective parametrization 
of power corrections that captures the dependence on NA/m for iV ^> m/A, where the 
power expansion completely breaks down. Even if one assumes that non-perturbative 
corrections are directly proportional to the corresponding renormalon residues, it remains 
difficult to control them in practice since the residues vary going beyond the large-/?o limit. 
Indeed, little is known about the Borel transform B[ag ucl ](u) in (16.141) beyond the large-/?o 
limit and away from the vicinity of the origin. When using the dispersive formulation 
instead, the problem translates into the parametrization of the Sudakov effective charge 
a| ucl (/x 2 ) itself in the infrared region. This should be an easier function to constrain, espe- 
cially if a| ucl (// 2 ) tends to a finite limit in the infrared already within perturbation theory, 
making it a slowly varying function. Constraining the coupling over the infrared region 
may be further helped by the universality of its infrared limit. It should be nevertheless 
stressed that the dispersive approach is advantageous only upon assuming that the rel- 
evant couplings are causal, free of Landau singularities. In the opposite case using the 
Borel formulation [54] is more straightforward. 

Finally, parametrization of the Sudakov effective charges over the infrared region pro- 
vides an alternative to the shape-function approach, one that matches smoothly onto the 
perturbative description, does not require the introduction of an explicit infrared cutoff, 
and most importantly, is better constrained. 
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A Coefficients of the Sudakov anomalous dimensions 



in MS 

The coefficients of the cusp anomalous dimension, corresponding to (12.91) are known to 
three-loop order [77]: 
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The coefficients of the jet-function anomalous dimension B (associated with an unre- 
solved jet with a constrained mass) are known to three-loop order [7]: 
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The coefficients of the Drell-Yan anomalous dimension T> UY (associated with large- 
angle soft radiation from two lightlike partons that annihilate to produce a non-colored 
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heavy object) are known to three-loop order [8,9]: 
d° Y = 



14 1 2 \C A 



8 7 ' 



116 inn 2 10. 1 

+ 109vr 2 + — C 3 + — 

81 3 S3 2p 



23 



180 



-7T 



517 2 245 , 2345\ 

7T 2 + Cs 

324 18 216 / 



1711 19, 

1 C3 H 

144 3 s 30 



vr 4 7T 2 



+ 



2/V 



223 2 11 4 



7T H 7T 

1296 360 



485 
1296 

11 



11 , 18821 209 . 

— 71 H (3 

24 1728 36 360 



7T CfCa 



(A.3) 



The coefficients of the anomalous dimension T> QD corresponding to the heavy-quark 
distribution function [86] as well as the heavy-quark fragmentation function are known to 
two-loop order [32]: 
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B Renormalon sum: Borel and dispersive represen- 
tations 



There exist two convenient representations of a dressed gluon which result in two different 
formulations of resummation formulae: 

• The Borel method: The running coupling, which include the effect of dressing, 
can be written as 



a: . 
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duT(u) 



Jo 



A 2 
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where is the gluon momentum and n(/c 2 ) is the vacuum polarization function, 
renormalized at fi 2 . The coupling is defined in the spacelike region — k 2 > and 
then analytically continued to the complex momentum plane. The superscript V 
stands for the V scheme (defined by the potential between two heavy quarks) where 
Ay = A 2 es n (cf. f)4.20p ). For one-loop running coupling, T(u) = 1. In this case one 
recovers 
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An all-order resummation of running-coupling corrections in a given quantity R{Q 2 ) 
with a single dressed gluon can thus be achieved by performing the momentum 
integration with the modified propagator 

1 _ 1 .RON 

-k 2 -io (-k 2 -ioy+ u ' y ' ' 

This procedure directly yields the Borel representation of the perturbative sum in 
the large-/?o limit in the form: 

oo 

W 2 )| largc/3o = jj duT{u) {Q 2 / A 2 )~ u B{u) . (B.4) 
o 

The dispersive method: The dispersive representation of the dressed gluon [45, 
49,61,62] takes the form 



aY{— k 2 ) a s (p 2 ) 1 1 J f p v (m 2 )dm 




7i 7i l + n(fc 2 ) A) I J m 2 -k 2 l + k 2 /A 2 v 



(B.5) 



where py (py < 0) is the spectral density function, defined by the discontinuity of 
the coupling on the timelike axis, 

/ 2\ Pot r v ( 2 . n w i l(3 a s (p 2 ) Im{n(m 2 )} 
py[wf) = — Im {< (-m - z0)/tt} = . (B.6) 

7T 7T 7T |1 + ll(m^J| 

Taking the time-like discontinuity of the one-loop running coupling (IB. 2ft one ob- 
tains: 

Mm ) I one-loop = l n 2( m 2/ A ^ + ?r 2 ' ( B " 7 ) 

The dispersive integral (the first term in the curly brackets in (IB.5j) ) then yields^ 

oo 

at PT (-k 2 ) 1 f dm 2 1 1 



7T poJ m?-k 2 \n 2 (m 2 /A 2 v ) + 7r 2 /3 \\n(-k 2 /A 2 v ) l + k 2 /A 2 v 

(B.8) 

that differs from the original one-loop coupling (1B.2j) by pure power terms that 
eliminate the Landau singularity. The second term in the curly brackets in ( IB. 51) 
cancels this additional term, and thus restores the Landau pole, making (1B.5I) con- 
sistent with ( IB. 21) . Note that we shall be using (1B.5|) rather than (1B.8|) : we will 
not assume anything about the way analytic properties of physical quantities are 
eventually restored in (non-perturbative) QCD. 

All-order resummation can be achieved using (IB. 51) by performing the momentum 
integration with a massive gluon propagator: 

1 1 



-k 2 — iO m 2 — k 2 — iO 



(B.9) 



^The superscript APT on (|B.8j) stands for 'Analytic Perturbation Theory' [88]; it distinguishes this 
object from the ordinary one loop coupling. 
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This results in the dispersive representation of the perturbative sum in the large-/3o 
limit: 

R (Q 2 )L Re0n = wA / ^(- 2 ) [HmVQ 2 )-m] + lH-^vlQ 2 ) - ^(0)] 

-a^{m 2 )T{m 2 /Q 2 ) + [f(-A 2 v /Q 2 ) - .F(O)] } , 

(B.10) 

where in the second line we applied integration-by-parts using (13.221) and 

T{m 2 IQ 2 ) = - m 2 A_T{m 2 IQ 2 ). 
dm z 

C Taking the large— N limits of finite— N characteristic 
functions 

In Sec. E]we identified the Sudakov limit of the characteristic functions in momentum space. 
A similar identification can be done in moment space. Let us now show that the Sudakov 
characteristic function in Eq. (I5.39j) can be recovered as the appropriate N —> oo limit of 
the moment-space characteristic function of the corresponding physical evolution kernel. 
We will also show that the formal 'jet' and 'soft' exponents (15.421) and (I5.43p . respectively, 
naturally emerge as limits of these kernels. Again we illustrate these statements in the 
examples of deep inelastic structure functions and the Drell-Yan cross section. 

Deep inelastic structure functions 

The finite-iV moment-space characteristic functionll 

£dis(f 2 ) (e, N) = / dxx N ~ 1 J r mS (F 2 )(e,x) (C.l) 
Jo 

appears in the dispersive representation (valid in the large-/?o limit) of the moment-space 
evolution kernel K ms(F2) (N,Q 2 ) (Eq. (ETBj) ): 

K ms{F2) (N, Q 2 ) =-f ^py(/i 2 ) [<3 ms{F ^ 2 IQ\ N) - g ms{F2) (0, N) 

large ft ft io I 1 V 

(C.2) 

This is the moment-space version of Eq. (I3.35p . Let us now split 9 mS (F 2 ) (e, N) into its real 
and virtual components: 



Q ms(F2) (e, N) = Cs(f 2 ) N) + G ( Z(f 2) W • (C3) 



II In Ref. [49] Qms(F 2 ) ( £ j N) is denoted by Tn (e), whereas ^dis(f 2 ) ( e ) an< ^ ^dy ( e ) are denoted by V s (e) 
and Vt (e), respectively. Note also that our normalization of the characteristic functions is half the one in 
this reference. 
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Considering first the real contribution, one finds [70], using the explicit expression in [49] 
for characteristic function: 

^ { ;l (F2) (e,N) = ^e- N ^j(e/r)^gj(Ne) , 
N -> oo Jo r (C.4) 

Ne fixed 

which is the analogue of the momentum space relation, Eq. (13.451) . Taking one derivative, 
we deduce: 



^€1 {F2) N)= - e~ Nr Tj (e/r) = Qj (Ne) 
N -> on Jo r 



N -> oo 

ATe fixed (C.5) 

1+ £_c^)k^« ) + (£_»)<-~. 

Let us turn now to the virtual corrections. Considering the same limit the virtual pieces 
Q( v > (e) as well as its derivative, (e) diverge. Indeed for e — > one gets, using the 
explicit expressions in [49]: 

^U( e )-4 ln2(e) -^ ln(e) -y-i (a6) 

and thus 

hm^ ( e )^ -In (iV) + In (iV6) + | + 0(1/N) . (C.7) 
Ne fixed 

Next, consider the combination £ D is(F 2 )( e ; N) — £ D is(f 2 )(0, N) occurring in the dispersive 
representation Eq. (1C.2|) and decompose it as 



g ms{F2) (e, N) - ^ DIS(F2) (0, N) = g£l {F2) (e, N) + (g^ m (e) - S dis (f 2 )(0, N)) . (C.8) 
Using now Eq. (4.43) in [49], one gets 

lim & ja(Fa) (p,N)~-hiN + ~-'y E , (C.9) 

and thus, using Eq. (IC.7I) : 

lim fe£> (e) - g ms{F2) (0, N)) = In (Ne) + \ + lE . (CIO) 
Ne fixed 
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Using Eq. (1C.5I) we deduce: 



lim £ DIS(F2) (e, N) - G ms{F2) (Q,N))=Qj (Ne) + In {Ne) + -t + 1e 
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and, comparing with Eq. (15.391) . we find: 



lim (g ms(F2) (e, N) - g ms(F2) (0, N)) = AQ { j +v) (Ne) . 
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Thus, as announced, the 'jet' Sudakov characteristic function, AQj +v ^ (Ne), is identified 
as the large- N limit with fixed iVe of the full finite- iV characteristic function of F 2 . Note 
also that, taking the derivative of Eq. (10. Ill) . one gets [68] 



limg ms(F2) (e,N) =Gj(Ne) - 1 
N -> oo 

Ne fixed 



(C.13) 



where —1 is the virtual contribution. Finally, taking the limit N —>■ oo inside the integrand 
of Eq. (lC.2j) (which yields an ultraviolet divergent integral) , and comparing with Eq. (15.451) , 
one obtains [70]: 
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Drell-Yan 

In a similar way the finite- A" characteristic function in the Drell-Yan case 
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appears in the moment-space equivalent of (13.391) : 
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Splitting the characteristic function into the real and virtual contributions we have: 

Guy (e, AT) = g$ (e, N) + gj$ (e) . (C.17) 
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For the real emission part, the large-A" limit with fixed N 2 e exists, and yields: 

lim<?£ (e, N) = 2 f° -e~ Nr ^ Y (e/r 2 ) = 2fo DY (iV 2 e) , (C.18) 
N -> oo -/o r 

N 2 e fixed 

which is the analogue of the momentum space relation, Eq. (13.461) . Taking a derivative we 
get: 

Km^ (e, N) = 2 / - e"^ ^ 5dy (e/r 2 ) = 2 fo DY (iV 2 e) = 2K (2^) . (C.19) 
N -> oo -/o r 

N 2 e fixed 

The virtual contribution to the characteristic function ^ Y (e) as well as its derivative, 
Q^"> (e) diverge in this limit. The small-e expansion is [49]: 

e ) ( e )-4 ln2(e) -^ ln(e) + y-i' (a20) 



so 



lim 3$ (e) ~ -2 In (N) + In (A^ 2 e) + - + 0(1/ N) . (C.21) 
N -> oo 2 
N 2 e fixed 

Decomposing the difference entering flC.161) according to 

£ DY (e, iV) - £ DY (0, iV) = 0&? (e, N) + (e) - £ DY (0, iV)) , (C.22) 
and using the relation 

(0, iV) = 2 

^DIS(F 2 ) 

(0,iV), (C.23) 

which follows from the results in Sec. 4.6 of [49] (and reflects the fact that the Drell-Yan 
cross section in the DIS scheme is an infrared and collinear safe quantity) one gets: 

lim £ DY (0, N) ~ -2 In iV + - - 2 lE + 0(1/N) . (C.24) 

N^oo 2 

Using Eq. (1C.21I) . one thus finds 

lim fe Y } (e) - £ DY (0, iV)) = In (iV 2 e) + 2 lE , (C.25) 
A^ oo y 
fixed 

which implies 

lim (g DY (e, N) - £ DY (0, N)) = 2 (q Sdy (N 2 e) + J In (iV 2 e) + Tb ) 
iV X oo y V ^ / 

N 2 e fixed (C.26) 

= 2 [K (2Nyft) + In (A^) + lE 
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and, comparing with Eq. (I5.39p . we find 

lim fe DY (e, N) - £ DY (0, N) 

N 2 e fixed 



(C.27) 



Thus, as announced the large-iV limit of the full finite-iV characteristic function with 
fixed Ne 2 reproduces the Sudakov characteristic function AGg^ (N 2 e). Note that in 
accordance with the discussion following Eq. f 1 5 . 2 2 j) above, the explicit result in Eq. (1C26|) 
is identical to the function appearing in Eq. (54) in [17] at b = 0. Note also that upon 
taking; one derivativdB of Eq. (10261) one obtains [68] 

lim£ DY (e,A0 = 2Qs m (N 2 e) - 1 . 
N -> oo 
N 2 e fixed 



(C.28) 



where —1 is the virtual contribution. Finally, taking the large- N limit inside the integral in 
Eq. (1C.16I) (which is yields an ultraviolet divergent integral) and comparing with Eq. (15.451) 
above one finds [70]: 



lim K UY (N, Q 2 
N -> oo 

Q 2 /N 2 fixed 



large f3 



E Sdy (Q 2 /N 2 ) 



large /3 



(C.29) 



D Renormalons and Landau singularities: example 



Let us consider for example the power terms distinguishing between ( 15. 9ft and ( 16. 14ft in 
the case of inclusive B decays under the assumption that the effective charges are given 
by (I4.39p . namely they do have Landau singularities. We obtain: 



In (^Sud(m 2 , N) 

_ C F r 2/m2 de 
A> Jo e 



Dispersive 



In Sud(m 2 ,iV) 



Borel 



Gj (eN) - Gj (e) ) - [Gs (eiV 2 ) - Q s (e) 



2(3 



j(N 2 - 1) -JVln(JV) - (N - 1) Hn 



soft 



jet 



A 2 

<:,: 1 +1E-5 



,,2 — _ A2 
r — yi GB 



A; 



77T 



O 



A 2 

GB 

m 2 



+ 



inCp 
2AT 



(N 



A 2 



1/2 



-(N-1) 



A 2 
m 2 



18 V 



soft 



jet 



soft 



A 



GB 



3/2 







A 2 

GB 

m 2 



**Note that Gs DY (N 2 e) does not satisfy the analogue of the Laplace representation Eq. (|5.3|) . owing 
the singular behavior of the momentum space characteristic function in the Drell-Yan case. 



1) 
to 



75 



where the first few terms in the expansion have been computed explicitly using (I5.13P 
and (15.171) . The first line summarizes the leading real power term that contributes to 
the (unambiguous) difference between ( 15.91) and the Principal Value of the Borel sum 
( 16. 14ft ; these terms are not related with renormalons. The second line summarizes the 
imaginary power terms which represent the ambiguity of the Borel sum owing to infrared 
renormalons. They provide an indication of the potential size of genuine non-perturbative 
effects. Note that these terms originate in the non- analytic terms in the small fi 2 expansion 
of the characteristic function. Obviously, there is one-to-one correspondence between 
these non-analyticities and the Borel singularities in ( I6.14p . The latter have been discussed 
in [40-43,54]. 

It is straightforward to identify the origin of the different terms in (ID.lj) . as indicated 
below each term: terms that scale at large iV as NA/m are associated uniquely with the 
'soft' (quark distribution) function, while those scaling as NA 2 /m 2 with the jet function. 
It is evident that both these classes of corrections contribute to both the real and the 
imaginary parts of (ID. ID . In practice, power corrections on the soft scale are important 
while those on the jet-mass scale can usually be neglected. Let us therefore consider the 
former in some detail: 

• The leading power term, 0(AN/m), represents the u — 1/2 renormalon ambiguity. 
This ambiguity has been shown [40] to cancel in the physical spectra with the am- 
biguity in defining the b quark pole mass, or equivalently A = Mb — in^, where Mb 
is the meson mass and m& is the b quark pole mass. 

In Refs. [40-43,54] the Principal Value Borel sum was used to define both the Su- 
dakov factor and the pole mass, eliminating the ambiguity. In Appendix [E] we provide 
a definition of the pole mass that suites the dispersive approach. 

• As a consequence of the Landau singularity in a^ ucl (k 2 ) assumed abova^l], a term of 
O ((AN/m) 2 ) appears as a real contribution to the difference between (I5.9P and the 
Principal Value of the Borel sum (16.141) . This is a parametrically large correction 
that directly influences the width of the computed spectra. It should be explicitly 
subtracted out when using the dispersive technique: it cannot be absorbed into the 
definition of any non-perturbative parameter, as there is no corresponding renor- 
malon ambiguity. 

• Finally, power terms of O ((AN/m) p ) where p is an integer p > 3, appear as renor- 
malon ambiguities. Corresponding non-perturbative corrections are expected to ap- 
pear as a consequence of the interaction between the b quark and the light degrees 
of freedom in the meson [40], distinguishing the quark distribution is the meson 
from that in an on-shell b quark. These power corrections have been parametrized 
in [41]. As far as these corrections are concerned the dispersive representation of the 
exponent (or its quark distribution function part) provides an alternative definition 
of the perturbative component, which is a priori as good as the definition based on 
the Principal Value Borel sum. 

^This correction does not appear if a| ucl (fc 2 ) has a causal analyticity structure, see [58] and the dis- 
cussion in Sec. [7] above. 
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To conclude, we have seen that 

• Genuine non-perturbative effects are related to renormalon singularities. As usual, 
the same conclusions with regards to such corrections can be reached using the dis- 
persive and the Borel formulations. As far as the regularization of the renormalons 
is concerned the two definitions are equivalent in principal: they differ by power cor- 
rections of the same parametric form one needs to introduce in order to parametrize 
genuine non-perturbative effects. Obviously, the non-perturbative parameters need 
to be defined according to the regularization used. 

• Quite independently of this, when the effective charges have Landau singularities 
the dispersive integral differs from the Borel sum by an additional set of computable 
power terms that are not related with renormalons; these include parametrically 
important contributions, 0(A 2 N 2 /m 2 ). 



E Definition of the pole mass in the dispersive ap- 
proach 

In Refs. [40-43, 54] the Principal Value Borel sum was used to define both the Sudakov 
factor and the pole mass, eliminating the ambiguity (see e.g. Eq. (4.4) in [43]). The 
use of the dispersion integral to define the Sudakov exponent requires a corresponding 
regularization of the pole mass (or A). Importantly, the relevant mass, which we shall 
refer to as the "dispersive pole mass", m^', is uniquely fixed by the formalism, and 
similarly to the Principal Value pole mass, it can be accurately determined based on some 
short distance mass (e.g. the MS mass). The relation between the "dispersive pole mass" 
and the Principal Value pole mass becomes transparent upon recalling the cancellation 
mechanism of the u — 1/2 renormalon. For concreteness let us refer to the B — > Xg'-f 
example. 

By combining the Laplace weight of the inverse-Mellin transform in (12.41) with the 
Sudakov factor (15.91) : 



2E. x 



7 \ 0,,J/™2 



Sud(m", N) 



x exp 




(E.l) 



where have used the relevant (square-root) term in the expansion of the characteristic 
function Gs Q r>{N 2 e) in (15.171) and omitted all other terms, which are irrelevant for the 
u = 1/2 renormalon. We consider the integration up to an arbitrary scale /i/ that should 
be above A; the specific scale choice will not affect the final result. 

The equivalent integral in the scheme-invariant Borel formulation can be obtained by 
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inserting (16.41) into (lE.lj) . which yields: 



2E*y 

m b 



N-l 



x exp 



Sud(m 2 ,iV) 
iVA PV fijNC F Tr 



u=l/2 



PV 



m b 



2E y 
Ms 

du 



N-l 



(E.2) 



l-2u \fj,j 



A 5 



2 \ u 



~~2 I T{u) 



sin iru 



B 



ITU 



l s QD 



where we have chosen the Principal Value prescription for both A and the u— 1/2 ambi- 
guity in the quark distribution function. 

The key point is that the object considered in ( IE. II) and ( IE. 21) is unambiguous, so the 
relation between the two definitions of the pole mass can simply be read off comparing 
the two equations: 



A 



disp. 



A r 



PV disp. 

rriu — m. 



-7716- 



du 



1 - 2u \fij 



A 2 



2 \ u 



dfj, 



2 -S^ 2 ) 




(E.3) 



-2) T(u 



sin TTU 



B 



nil 



Moreover, the r.h.s. of Eq. (IE. 31) does not depend on the cutoff fii (see Eq. (B.19) in [58] 
and Sec. 3.2 in [45]). We can thus set \ii — A to get: 



A 



disp. 



At 



.PV 



+ ^PV 




2 <i k (^ 2 : 




A 2 



(E.4) 



u 



exhibiting the fact that this difference is a number of order A which is independent of any 
hard scale such as fii or m b . Note that these equations are valid even if the Euclidean 
coupling has Landau singularities (such as the one loop coupling). The Minkowskian 
coupling is still infrared finite in this case, and the convergence of the Borel integral is 
ensured by the oscillations of the factor sin(7ru)/7r-u. 

It is also possible to give a Euclidean version of the previous definition, iff one assumes 
the Euclidean Sudakov effective charge is causal and infrared finite. It is straightforward 
to show that the analogues of Eqs. (IE. II) . (IE.3I) and (IE.4I) are: 



f2E. 



N-l 



Sud(m\N) 




l N 2 k 2 



mf 



(E.5) 
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a a pv disp. CV f^ 1 dk 2 B j 2 \ 

Adisp. - A PV = m b - m b = -m b — / — a SQD (k ) J 

° V b (E.6) 

C p f°° du fA 2 \ u m , . „ 



2* 



a Eucl 



and, setting \ij = A, since the r.h.s. of Eq. (IE. 61) does not depend on the cutoff \ij in the 
case of a causal coupling (see Eq. (3.19) in [58]): 



A dls , - A PV = mr - mf- = A x j - 9l J* ^ a - (fc2) ^ 



(E.7) 



Note also the convergence of the Borel integral has to be insured by oscillations in 
B 



-.Eucl 
^QD 



m) if CLsqn(k ) has an infrared fixed point (similarly to the Minkowskian cou- 
pling above). Finally, we note that, although Eqs. (IE. II) and (IE. 51) define different u— 1/2 
pieces of the Sudakov exponent, the corresponding mj lsp ' masses are the same: this result 
follows immediately comparing Eq. (3.19) and (B.19) in Ref. [58]. 
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